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ABSTRACT 


This  report  presents  the  results  of  a  study  made  by 
Monsanto  Research  Corporation  during  the  final  phase  of 
Contract  No.  DAAH01-72-C-0490.  The  intent  of  the  study  was 
to  define  those  physical  properries  of  magnetic  bubble 
materials  which  need  to  be  characterized,  to  investigate 
the  known  methods  for  making  these  characterizations,  and, 
finally,  to  recommend  those  characterization  procedures 
which  appear  "best"  at  this  time.  Those  material  properties 
included  in  the  study  were:  Substrate  and  film  composition  - 
substrate  and  film  lattice  parameter  -  substrate  and  film 
defect  identification  and  location  -  film  thickness  -  film 
thickness  variations  -  characteristic  length  and  domain 
dimensions  -  saturation  magnetization  and  magnetic  fields  - 
domain  wall  energy  -  domain  wall  dynamic  properties  - 
anisotropy  -  coercivity  -  magnetostriction  coefficients  - 
reorientation  temperature  -  Neel  temperature  -  compensation 
temperature  -  temperature  variation  of  material  parameters  - 
thermal  conductivity  -  refractive  index  -  Faraday  effect  - 
Kerr  effect  -  optical  and  magneto-optical  absorption  - 
electrical  conductivity  -  photo  conductivity. 

Extensive  reference  to  existing  literature  has  been 
made  in  an  endeavor  to  keep  this  report  within  reasonable 
size.  As  a  result,  a  lengthy  bibliography  has  been  compiled 
to  which  the  reader  is  urged  to  refer  for  a  more  detailed 
analysis  of  the  measurement  procedures  discussed  here. 
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1.  INTRODUCTION 

The  concept  of  a  magnetic  memory  operating  without  moving 

parts  has  been  of  interest  to  the  computer  industry  for  many 

years.  Indeed,  the  ferrite  core  memories  in  wide  use  today 

are  such  devices  but  the  complexity  and  resulting  cost  of 

their  assembly  has  kept  interest  in  other  alternatives  high. 

A  number  of  memory  related  concepts  involving  the  motion  of 

magnetic  domains  in  metallic  ferromagnetic  materials  have 

reached  the  working  device  level  in  the  form  of  shift 

TM 

registers,  including  Hughes  Aircraft  Dynabit  and  Cambridge 

TM  • 

Memories,  Inc.  DOT  Systems.  Bell  Telephone  Laboratories 

has  most  strongly  pursued  the  concept  of  moving  domain 

memory  devices  in  insulating,  ferrimagnetic  materials, 

starting  first  with  the  orthoferrites  (prototype  -  YFe03) 

and  progressing  more  recently  to  the  magnetic  garnets 

(prototype  -  Y3Fe5°i2^  *  In  these  materials  the  desirable 

configuration  is  one  in  which  the  unique  easy  axis  of 

magnetization  is  perpendicular  to  a  thin  slab  or  film  of  the 

material.  In  the  presence  of  a  magnetic  bias  field  of 

the  proper  magnitude  along  the  magnetization  direction, 

cylindrical  magnetic  domains  of  opposite  magnetic  polarity 

will  be  stable  in  these  materials.  These  cylindrical  domains, 

when  viewed  between  crossed  polarizers,  look  like  tiny 

bubbles  by  virtue  of  the  Faraday  effect  in  the  magnetic 

material.  This  term  "magnetic  bubble"  has  become  widely 

used  to  describe  a  cylindrical  magnetic  domain.  A  great  deal 

of  work  has  already  been  done  but  much  work  remains  in  order 
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to  be  able  to  generate,  move,  and  detect  bubbles  with  the 
speed  and  reliability  required  by  the  computer  industry. 

Many  people  are  now  convinced,  however,  that  this  can  and  will 
be  done.  Still  further  development  will  be  required  to  bring 
to  fruition  the  potential  for  logical  operations  which  is 
inherent  in  the  repulsive  interaction  which  exists  between 
the  bubbles.  The  appeal  of  a  complete,  if  small,  bubble 
computer  is  considerable  and  such  efforts  will  undoubtedly 
continue. 

A  very  significant  part  of  the  development  of  magnetic 
bubble  technology  has  been  and  will  continue  to  be  the 
characterization  of  the  properties  of  appropriate  magnetic 
materials.  Older  measurement  methods  are  being  adapted  to 
the  special  requirements  of  these  thin,  transparent  magnetic 
bubble  films,  and,  in  addition,  new  methods  of  measurement 
are  being  developed  which  take  advantage  of  the  unique  pro¬ 
perties  of  these  films,  An  particular  the  easily  observed 
Faraday  effect  in  the  visible  portion  of  the  electromagnetic 
spectrum.  There  is  no  doubt  that  further  refinements  will 
take  place  in  the  methods  and  the  properties  used  to 
characterize  these  materials.  Improved  accuracy,  more 
reliable  prediction  of  device  performance,  and  more  rapid 
data  acquisition  and  analysis  are  all  desirable  and  wj 11  be 
achieved. 

This  contract  requires  the  conduct  of  a  study  of  the 
state  of  the  art  in  the  characterization  of  single-crystal, 
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magnetic  films, 


non-magnetic  substrates  and  single-crystal  thin 
including  establishment  of  all  the  items  required  to  be 
characterized,  presentation  of  the  known  techniques  for  each 
of  the  methods  of  characterization  and  identification  of  the 
recommended  characterization  method  for  each  item.  The 
present  report  has  been  prepared  in  fulfillment  of  that 
requirement.  We  have  relied  upon  references  to  the 
available  literature  wherever  possible  as  the  traditional 
and  best  means  of  presenting  much  of  the  detail  of  these 
techniques.  In  order  to  increase  the  usefulness  of  the 
report  we  have  also  included  references  to  the  application 
of  the  various  techniques  to  magnetic  bubble  domain 
materials  wherever  possible.  For  the  present  purposes 
these  materials  are  defined  as  the  orthoferrites  and  iron- 
containing  garnets,  with  the  latter  receiving  much  of  the 
attention  because  of  their  more  desirable  properties  — 
primarily  smaller  bubble  size  and  ease  of  epitaxial  growth. 

The  measurement  of  the  following  properties  of  these  materials 
are  discussed: 

Composition 

Lattice  constants 

Substrate  defects 

Film  thickness  and  its  variation 

Film  defects 

Characteristic  length 

Magnetization 


3 


Domain  wall  energy 
Domain  wall  dynamics 
Ferromagnetic  resonance 
Anisotropy 
Coercivi ty 

Magnetostriction  coefficients 
Thermal  properties 
Optical  properties 
Electrical  properties 

Not  all  of  these  items  are  felt  by  the  authors  to  be 
worthy  of  routine  measurement,  but  are  discussed  here  more 
for  reasons  of  completeness  than  for  any  reflection  on 
their  importance.  The  authors  of  this  report  are  Roger  W. 
Shaw,  Robert  M.  Sandfort,  and  Jerry  W.  Moody.  They  would 
welcome  comments,  suggestions,  and  criticism  from  the 
readers . 
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2.  COMPOSITION 


A  knowledge  of  the  chemical  composition  of  the  substrate 
as  well  as  the  epitaxial  layer  is  required  for  the  complete 
characterization  and  understanding  of  garnet  bubble 
materials.  Impurities  or  stoichiometric  deviations  in  the 
substrate  can  affect  the  lattice  constant  of  the  material. 

Even  small  variations  of  lattice  constant  can  result  in 
significant  effects  in  the  epitaxial  film.  This  is 
especially  true  for  films  exhibiting  stress -induced 
anisotropy  for  which  the  non-cubic  anisotropy  depends  upon 
the  strains  resulting  from  a  deliberate  mismatch  of  the 
lattice  constant  of  film  and  substrate.  In  addition,  the 
chemical  components  of  the  substrate  can  be  physically 
incorporated  into  films  grown  by  either  chemical  vapor 
deposition  or  by  liquid  phase  epitaxy.  For  example, 

Robinson,  et  al  (Ref.  1  ),  observed  an  interfacial  layer 
between  Tb2> 5ErQ >5Fe5012  films  grown  by  CVD  on  Sn^Ga^ 
substrates.  They  reasoned  that  the  intermediate  layer 
contained  some  samarium  from  the  substrate.  In  the  LPE 
growth  of  garnet  films  in  PbO-based  solvents,  there  is 
ample  evidence  that  the  surface  of  the  substrate  is 
dissolved  before  film  growth  is  initiated.  The  dissolved 
components  can  then  be  reprecipitated  as  the  film  grows. 

Levin ste in f  et  al  (Ref.  2)  observed  intermediate  layers  in  LPE 
films  of  Gd3_xTbxFe5012  on  Nd3Ga5012  substrates.  The 
components  of  the  solution  itself  are  another  source  of 


impurities  which  can  be  incorporated  in  the  film.  Since 
the  magnetic  properties  of  the  film  are  a  sensitive 
function  of  the  composition,  it  is  important  to  know  not 
only  the  major  constituents  but  also  what  contaminants 
might  be  present  in  the  film. 

2 . 1  Film  Composition 

The  problems  associated  with  the  determination  of  the 
composition  of  magnetic  garnet  films  have  been  discussed 
by  Mee,  et  al  (Ref.  3  ).  In  general,  the  films  of  interest 
are  less  than  lOym  thick  and  only  a  minute  quantity 
is  available  for  analysis.  Usually  both  the  film  and 
substrate  contain  at  least  gallium  in  common  and  separation 
of  film  and  substrate  is  not  a  trivial  matter.  In  typical 
multi- component  films  the  physical  and  magnetic  properties 
of  the  films,  such  as  lattice  constant,  magnetization, 
etc.,  are  not  specific  and,  thus,  can  not  be  used  as  a  basis 
for  a  reliable  determination  of  composition.  Most  of 
these  problems  are  avoided  by  electron  probe  microanalysis 
(EPM)  for  compositional  determinations  (Ref.  4  ) . 

EPM  offers  three  advantages  over  other  means  of 
compositional  analysis. 

1)  The  effective  volume  of  characteristic  X-ray 
production  is  small,  typically  of  the  order  of  2ymJ. 

This  shallow  penetration  of  the  epitaxial  layer  insures 

that  only  the  area  of  interest  is  considered  in  the  determination. 

2)  Due  to  the  small  area  of  analysis,  it  is  possible 
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to  perform  homogeniety  determinations  over  the  surface  of 
the  epitaxial  layer. 

3)  Since  the  technique  is  non-destructive  it  is 
possible  to  utilize  the  film  in  further  testing  and  evaluation 
subsequent  to  EPM  analysis. 

The  concentrations  of  all  elements  composing  the 
garnet  matrix  are  determined  quantitatively  with  the  excep¬ 
tion  of  oxygen.  The  concentration  of  oxygen  in  all 
determinations  is  assumed  to  compose  the  balance  of  the 
composition.  This  is  necessary  since  the  instrument 
currently  used  at  Monsanto  for  the  analysis  does  not  have 
the  capability  to  detect  elements  lighter  than  aluminum. 

For  the  samples  evaluated  to  date,  the  concentration 
of  rare  earth  metals,  lead,  and  iron  have  been  determined 
by  comparison  to  pure  metal  forms  of  these  elements.  The 
determination  of  gallium  concentration  was  made  by  comparison 
to  gallium  phosphide.  The  purity  of  the  standards  used 


is  as  follows: 

Rare  Earths 

-  99.9% 

Fe 

Spectrographic  grade 

Pb 

Spectrographic  grade 

GaP 

-  Monsanto  Semiconductor 

Matrix  corrections  are  applied  when  the  sample  counts 
are  compared  to  the  standards.  The  accuracy  of  the  method 
is  approximately  3%  of  the  actual  concentration  of  the  metals. 
Since  the  electron  beam  samples  an  area  only  about  2 
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microns  in  diameter,  the  compositional  uniformity  of  the 
epitaxial  layers  is  determined  by  taking  counts  at  a  number 
of  widely  separated  points  on  each  layer.  The  agreement 
of  these  determinations  is,  in  general,  well  within  the 
experimental  accuracy  of  the  method.  The  compositioral 
uniformity  of  the  layers  is  also  checked  by  measuring 
the  Neel  temperature  at  various  points.  This  measurement 
is  discussed  in  a  later  section.  Again,  the  results  at 
various  points  on  a  given  sample  generally  agree  within 
the  accuracy  of  the  method. 

2 . 2  Trace  Impurity  Analysis 

The  oxides  (rare  earth,  iron,  gallium,  lead,  etc.)  used 
in  garnet  substrate  and  film  growth  are  now  available  in 
high  purity  from  a  number  of  sources.  Also,  the  growth 
of  crystals  by  the  Czochralski  method  generally  effects 
a  further  purification  (Ref.  5  ).  Thus,  impurities  in 
garnet  substrates  are  usually  present  in  trace  amounts 
(with  the  exception  of  iridium  in  the  substrate  and  lead 
in  flux  grown  films)  and,  to  date,  have  not  been  considered 
a  serious  problem  in  garnet  bubble  material  technology. 

This  attitude  may  change  as  the  technology  advances  and 
ever  more  stringent  requirements  are  placed  on  compositional 
control. 

The  analysis  of  rare  earth  garnets  (or  starting 
materials)  for  impurities  is  relatively  straight  forward 
and  may  be  accomplished  by  a  number  of  well  established 
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methods.  These  include  emission  spectroscopy,  mass 
spectroscopy,  chemical  spectrophotometry  and  fluorometry, 
and  atomic  absorption  (see  references  6 
and  7  for  description  of  these  and  other 

methods  for  trace  analysis) .  It  should  be  pointed  out 
that  these  methods  are  destructive  but  usually  require  only 
a  small  amount  of  sample.  All  these  methods  are  capable 
of  accurate  determination  of  trace  amounts  and  the 
choice  of  a  particular  method  depends  primarily  on  the 
availability  of  specific  equipment  or  personnel. 

The  analysis  of  bubble  materials  for  trace  impurities 
has  been  seldom  reported.  However,  semiquantitative 
analysis  by  emission  spectroscopy  has  been  performed  at 
Monsanto  on  selected  crystals  and  materials  which  were 
believed  to  be  contaminated  by  impurities.  The  technique 
used  is  that  discussed  by  Wang,  et  al  (Ref.  8  ). 

2.3  Stoichiometric  Deviations 

The  compositional  range  of  stability  of  the  rare  earth 
garnets  is  very  narrow.  In  addition,  it  is  probable 
that  the  excess  components  are  not  distributed  uniformly 
but  are  concentrated  in  the  growth  striations  and  "core" 
region  of  Czochralski  grown  crystals.  It  is  difficult 
to  determine  stoichiometric  deviations  accurately  in  such 
systems  so  that  completely  satisfactory  analytical 
techniques  have  not  yet  been  developed. 
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The  electron  microprobe  has  been  used  to  detect 
variations  in  the  gallium  concentration  in  wafers  of 
Gd3Ga2°12‘  Recently,  Braginski ,  et  al  (Ref.  9  )  used  the 
elect  on  microprobe  to  study  the  compositional  variation 
in  Y3Feij012  films  grown  by  chemical  vapor  deposition  on 
Gd3Ga5°12  substrates.  However,  this  method  has  not  been 
made  quantitative. 

Brandle,  et  al  (Ref.  10)  have  related  the  lattice 
constant  of  Gd^Gaj-O^  to  the  stoichiometry  of  the  melt  from 
which  the  crystals  were  pulled.  With  this  information 
it  is  possible  to  obtain  a  qualitative  estimation  of  the 
stoichiometry  of  a  Gd^Ga^O^  crystal  by  a  precise  lattice 
constant  measurement. 
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3. 


LATTICE  CONSTANTS  AND  THERMAL  EXPANSION  COEFFICIENTS 


An  accurate  knowledge  of  the  lattice  constant  of  both 
film  and  substrate  is  one  of  the  most  important  factors  in 
the  successful  growth  of  epitaxial  layers  and  in  the  inter¬ 
pret?  -ion  of  the  magnetic  properties  of  the  films.  This  is 
especially  true  for  those  film  compositions  which  exhibit 
stress-induced  anisotropy.  Film  and  substrate  must  be  matched 
closely  enough  to  yield  smooth,  uncracked  layers  with  a 
high  quality  interface,  but  mismatched  enough  that  bilateral 
stresses,  acting  through  magnetostriction,  result  in  the 
desired  anisotropy.  The  lattice  constant  of  Czochralski- 
grown  garnets  depends  upon  the  conditions  of  growth  and  will 
be  found  to  vary  over  a  relatively  narrow  range.  However, 
the  range  of  f ilm/substrate  mismatch  that  can  be  tolerated 
is  so  small  that  the  lattice  constant  of  each  substrate 
crystal  must  be  measured  and  the  composition  of  the  film 
adjusted  correspondingly. 

The  general  methods  for  determining  the  structure  and 
lattice  constant  of  crystalline  substances  are  well  known 
and  need  little  elaboration  here.  The  reader  is  referred 
to  any  of  the  many  excellent  texts  on  the  subject  (for 
example,  ref.  11).  It  is  appropriate,  however,  to  outline 
the  specific  procedure  used  at  Monsanto  to  determine  the 
precise  unit  cell  dimensions  of  the  substrates  and  films, 
as  an  example. 


11 


The  lattice  constant,  aQ  is  determined  by  the  x-ray 
diffractometry  procedure  published  by  Pierron  and  McNeely 
(Ref.  12).  Binoe  the  original  precision  of  the  method,  + 
O.OOOOiA,  is  not  needed  in  this  study  corrections  for  the 
Lorentz-polarization ,  x-ray  refraction  and  temperature 
errors  are  not  made.  The  precision  of  the  procedure  used 
and  described  herein  is  +  0.0001A. 

The  (111)  oriented  garnet  crystal  is  mounted  on  the 
sample  holder  and  the  (888) ,  (12  12  12) ,  and  (16  16  16) 
reflections  are  slowly  scanned  using  Mo  radiation.  A 
scanning  speed  of  0.125W1,  corresponding  to  one  inch 
per  minute  of  recorder  chart  travel,  permits  accurate  estimate 
Of  the  29  angle  at  maximum  peak  intensity  to  within  0.002*. 

The  d  spacings  are  computed  from  the  29  angles  and  are  then 
converted  to  corresponding  a„  values  for  the  (888) ,  (12  12  12) 
and^,l6  16  16)  reflections.  The  aQ  values  are  plotted  versus 
cot  9  so  that  an  extrapolation  to  high  9  angles  may  be  made 
which  will  minimise  systematic  misalignment  error.  The 
extrapolated  value  of  a „  at  cot29  -  0  yields  the  cell 
parameter,  if  the  crystal  has  been  aligned  properly,  the 
plot  yields  a  straight  line  parallel  to  the  cot29  axis,  m 
routine  practice  the  alignment  is  considered  satisfactory  if 
the  difference  in  a0  values  between  the  (888)  and  (16  16  16) 
reflections  is  no  greater  than  0.001A.  Alignment  is  routinely 
verified  by  slowly  scanning  the  (888)  and  (12  12  12) 
reflections  of  a  standard  GaAs  crystal,  aQ  =  5.65399.  The 
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2©  values  of  peak  maximum  should  be  60.23  and  97.63  degrees 
respectively  for  Mo  Ka^  radiation.  Fig.  1  is  a  typical  plot 
for  a  GGG  crystal.  In  this  case  the  curve  extrapolates  to 


an  a  value  of  12.3801A. 
o 

The  same  procedure  is  used  to  determine,  in  situ,  the 
lattice  constant  of  the  epitaxial  films.  In  this  case  two 
peaks  are  obtained:  one  peak  corresponding  to  the  strained 
lattice  constant  of  the  film  and  a  second,  usually  larger  in 
amplitude,  corresponding  to  the  strained  lattice  constant  of 
the  substrate.  The  lattice  constants  obtained  in  this  manner 
correspond  to  interplanar  spacings  which  are  perpendicular  to 
the  film  surface,  and  in  general  are  different  from  the 
unstrained  lattice  parameters.  The  lattice  dimensions  parallel 
to  the  film  surface,  which  are  of  more  interest,  can  be 
obtained,  according  to  Mee,  et  al  (Ref.  13)  by  the  relation¬ 
ship  between  parallel  and  perpendicular  strain/ 


_  (1  +  H)  _  v 

Aa  (1  -  v0  (as  af) 


(3-1) 

where  Aa  is  the  difference  in  lattice  dimensions  obtained 
from  the  x-ray  analysis,  y  is  Poisson's  ratio  for  the  film 
composition  and  ag  and  a^  are  the  unstrained  lattice  constants 
of  substrate  and  film,  respectively. 

The  authors  feel  there  is  some  question  as  to  whether 
one  actually  obtains  the  unstrained  lattice  constant  of  the 
film,  af,  from  this  analysis.  In  the  data  obtained  as  outlined 
above,  the  substrate  lattice  constant  determined  in  the 
presence  of  a  thin  epitaxial  film  is  different  in  general  from 
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its  value  determined  separately.  It  is  usually  shifted  in 
the  direction  of  the  film  lattice  parameter  in  agreement 
with  Shick's  (Ref.  14)  observations.  This  may  indicate 
that  a  kind  of  "average"  lattice  parameter  is  being  measured 
for  both  film  and  substrate,  or  that  the  substrate  is  in 
fact  strained  itself  near  the  interface.  We  feel  that  Aa'* 
is  an  important  parameter,  but  are  disinclined  to  interpret 
an  unstrained  film  lattice  parameter  from  Aa^  on  the  basis 
of  results  seen  to  date.  Probably  fundamental  lattice 
parameters  of  this  type  should  be  taken  from  bulk  single 
crystals  of  identical  composition,  or  from  very  thick 
epitaxial  layers  which  are  strong  enough  to  be  free 
standing. 

For  an  apparatus  equipped  with  a  high  temperature  cell, 
the  measurement  of  precise  lattice  constants  at  elevated 
temperatures  provides  a  convenient  means  of  determining  the 
thermal  expansion  coefficients  of  the  garnets  (Ref.  15) .  If 
such  equipment  is  not  available,  the  more  conventional 
dilatometry  may  be  used  (Ref.  16).  This  is  the  method 
used  by  Mee,  et  al  (Ref.  17). 


4 .  SUBSTRATE  DEFECTS 


Although  crystals  of  Gd^Ga,.  °12  can  now  be  grown 
essentially  free  of  all  defects,  the  need  for  substrate 
examination  and  evaluation  has  not  been  obviated.  The 
growth  of  epitaxial  magnetic  garnet  films  of  device  quality 
requires  substrates  of  the  very  highest  quality.  Indeed, 
most  of  the  defects  which  impede  domain  motion  in  the 
magnetic  film  can  be  traced  to  defects  already  present  in 
the  substrate  before  film  deposition.  Therefore,  substrate 
evaluation  will  always  remain  an  important  part  of  bubble 
material  technology. 

Substrate  defects  can  be  introduced  either  during  the 
crystal  growth  process  or  in  subsequent  cutting,  lapping 
or  polishing.  Scratches,  dust  and  solvent  stains  are 
examples  of  defects  which  fall  into  the  latter  class. 
Frequently  they  can  be  detected  by  close  visual  examination 
of  the  substrates  and  can  be  avoided  by  good  housekeeping 
practices  and  proper  polishing  procedures.  They  need  not 
cause  the  complete  rejection  of  the  substrate  since  the 
wafer  can  be  repolished  or  recleaned.  Defects  which  are 
introduced  during  crystal  growth  are  more  serious.  These 
include  strains,  inclusions,  voids,  growth  striations  and 
dislocations.  Experience  indicates  that  defect-free 
epitaxial  films  cannot  be  grown  on  substrates  containing 
such  defects.  Recently,  Chaudhari  (Ref.  18)  discussed  the 
defects  in  garnet  substrates  that  influence  magnetic  domain 
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motion.  The  detection  and  evaluation  of  these  defects 
are  discussed  in  the  following  paragraphs.  It  will  be 
apparent  that  most  of  the  techniques  recommended  for 
substrate  evaluation  depend  upon  the  interaction  of  defects 
with  light  or  X-rays,  Although  the  techniques  are  discussed 
in  reference  to  a  specific  defect  for  which  they  are  especially 
appropriate,  their  application  can  be  more  general.  For  example, 
X-ray  topography  can  not  only  be  used  to  study  striations 
and  banding,  but  it  also  reveals  dislocations  and  inclusions. 
Other  examples  of  complementary  methods  and  how  they  can 
be  used  to  advantage  will  be  obvious. 

4 . 1  Strain 

The  rare-earth  gallium  garnets  are  grown  by  the 
Czochralski  method  from  RF  heated  iridium  crucibles  at 
approximately  1800°C.  Large  thermal  gradients,  both  in 
and  above  the  melt,  are  natural  in  such  systems.  As  a 
result  of  the  thermal  stresses,  the  garnet  crystals  are  often 
strained.  Indeed,  one  of  the  most  characteristic  features 
of  Czochralski  grown  rare-earth  garnets  is  a  highly 
strained,  central  "core"  which  is  associated  with  the 
development  of  facets  at  the  solid/liquid  interface  during 
crystal  growth.  (The  core  can  be  eliminated  by  adjusting 
the  growth  conditions.) 

Strain  induces  optical  heterogeneity  and  bire¬ 
fringence  in  the  crystal  which  can  be  detected  by  examination 
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with  polarized  light.  When  the  crystal  is  viewed  between 
crossed-polarizers ,  the  strain-induced  birefringence  manifests 
itself  as  color  and  intensity  variations.  Fig.  2  illustrates 
the  appearance  of  a  strained  Gd^GagO^  crystal  when  viewed 
between  crossed  polarizers.  The  presence  of  the  central 
core  is  readily  apparent.  For  this  examination,  the  ends 
of  th  i  boule  must  be  cut  and  polished  parallel.  This 
qualitative  examination  for  strain  is  used  routinely  at 
Monsanto. 

A  more  quantitative  evaluation  of  the  strain  in  garnet 
crystals  can  be  obtained  by  Twyman-Green  interferometry 
(Ref.  19).  This  method  has  been  widely  used  for  assessing 
strain  in  sapphire,  ruby  (Ref.  20),  garnet  (Ref.  21)  and 
other  transparent  crystals.  The  principal  of  the  Twyman- 
Green  interferometer  is  shown  in  Fig.  3.  The  test  crystal 
is  mounted  on  one  arm  of  the  apparatus  and  the  interference 
pattern  is  recorded  by  a  camera  located  on  another  arm.  A 
variation  of  one  fringe  on  the  interferogram  of  a  crystal 
of  length  L  with  plane  parallel  end  faces  is  due  to  an 
optical  path  difference  of: 

A(nL)  =  -j-  (4-1) 

where  n  is  the  refractive  index  and  X  is  the  wavelength  of 
the  light  used. 

It  should  be  pointed  out  that  the  usual  Twyman-Green 
interferogram  integrates  the  effect  of  optical  inhomogeneities 


Picture  2. 
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crysLil  between  crossed  polarizers. 
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Figure  3.  Schematic  diagram  of  Twyman-Green 

interferometer.  After  Nestor,  ref.  23. 


along  the  entire  length  of  the  rod  so  that  the  specific 
source  or  location  of  the  inhomogeneity  cannot  be  identified. 
However,  more  specific  information  can  be  obtained  by 
immersing  the  crystal  in  an  index  matching  liquid  (contained 
in  a  cell  with  flat  parallel  faces)  and  viewing  from  the  side. 

4.2  Inclusions 

The  rare-earth  gallium  garnets  are  not  stable  at  their 
melting  points  in  an  inert  atmosphere.  The  melts  tend  to 
loose  gallium  via  a  volatile  suboxide  (Ref.  22).  The 
decomposition  of  Ga^  can  be  suppressed  by  maintaining  a 
partial  pressure  of  oxygen  above  the  melt.  However,  the 
garnet  melts  are  usually  contained  in  iridium  crucibles 
which  are  readily  oxidized  when  in  an  oxygen  atmosphere 
at  elevated  temperatures.  It  is  difficult  to  maintain 
melt  stoichiometry  and  minimize  crucible  oxidation 
simultaneously.  During  the  long  periods  of  time  required 
to  grow  the  crystals,  the  melt  tends  to  become  non- 
stoichiometric  and  saturated  with  Ir02-  It  is  not 
surprising,  therefore,  that  inclusions  of  one  or  more 
foreign  phases  are  a  common  defect  in  garnet  crystals. 

Brandle ,  et  al  (Ref.  10)  have  discussed  the  origin  and 
appearance  of  the  most  common  inclusions  found  in  Gd3Ga5012 
crystals.  Three  types  of  inclusions  were  identified: 

(1)  triangular  or  hexagonal  platelets  of  metallic  iridium 

(2)  dark  square  or  cubic  crystals  which  are  believed  to 
be  particles  of  a  gadolinium  gallium  suboxide, 
probably  gadolinium  orthogalliate  (GdGaOy 
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(3)  transparent  acicular  crystals,  tentatively  identified 
as  Gd2C>3. 

The  first  two  types  are  by  far  the  most  predominant,  with 
the  third  found  only  occasionally. 

Inclusions  in  the  garnets  can  be  studied  most  conveniently 
by  optical  microscopy.  Slices  of  the  crystal  polished  flat 
and  parallel  are  required  for  these  studies.  Examples  of 
the  most  common  defects  are  illustrated  in  the  photo¬ 
micrographs  of  Fig.  4.  The  density  of  the  inclusions  can 
be  obtained  by  counting  the  number  of  inclusions  in  a 
known  volume  of  the  wafer. 

Tyndall  scattering  may  also  be  used  to  measure  the 
density  of  inclusions  (Ref.  20).  Here  collimated  monochromatic 
light  is  directed  through  the  crystal,  usually  along  the 
growth  or  optical  axis.  The  emergent  scattered  light  (low 
angle  scattering)  can  be  measured  by  placing  a  film  strip 
in  an  arc  whose  center  is  the  center  of  the  exit  end. 
Alternatively,  the  large  angle  scattering  can  be  measured 
by  photographing  the  light  scattered  perpendicular  to 
the  collimated  beam.  In  this  mode,  the  crystal  can  be  scanned 
to  determine  the  site  and  distribution  of  scattering 
centers.  Nestor  (Ref.  23)  has  described  an  apparatus  to 
measure  large  angle  scattering  for  ruby  laser  rods . 

4 . 3  Dislocations 

The  energy  to  create  dislocations  in  the  rare  earth 
garnets  is  large  and  crystals  can  be  grown  essentially  free 
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Figure  4a.  Iridium  inclusions  (200X) 


Figure  4b,  Gadolinium  gallium  suboxide  particle  (1000X) 
Figure  4.  Common  inclusions  in  Gd3Gd5012 
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of  all  dislocations.  However,  dislocations  are  not  uncommon. 

Dislocations  will  affect  the  optical  homogeneity  of  crystals 

and  their  presence  can  be  inferred  by  examination  in 

crossed  polarizers.  For  example.  Fig.  5  is  a  photograph 

4  2 

taken  of  a  crystal  containing  about  10  dislocations/cm 

between  crossed  polarizers.  This  photograph  may  be  compared 

with  that  of  Fig.  2  which  is  of  a  crystal  containing 

2 

less  than  10  dislocations/cm  .  The  strain  in  the  low 
dislocation  density  crystal  is  concentrated  in  the  core 
region  whereas  it  appears  to  be  distributed  throughout  the 
high  dislocation  density  crystal. 

Dislocations  can  be  revealed  by  selective  etching. 
Orthophosphoric  acid  at  160-170°C  is  a  suitable  etchant 
for  the  rare  earth  gallium  garnets  (Ref.  24).  On  the  {111} 
planes  (the  most  commonly  used  plane  for  epitaxial  growth) , 
the  dislocation  etch  pits  are  six  sided  with  overall  triangular 
symmetry.  A  typical  dislocation  etch  pit  is  shown  in 
Fig.  6.  This  photomicrograph  was  made  with  the  optical 
microscope;  however,  the  morphology  of  the  pits  would  be 
more  apparent  under  the  scanning  electron  microscope. 

4 . 4  Growth  Striations 

Czochralski  grown  garnet  crystals  usually  exhibit 
distinctive  growth  striations.  It  is  believed  that  the 
striations  are  slight  compositional  variations  (stoichiometric 
deviations  and/or  impurities)  resulting  from  thermal 
fluctuations  in  the  melt  during  growth.  If  such  is  the 
case,  then  corresponding  variations  in  lattice  constant 
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Figure  5.  Gd3Ga5012  crystal  of  high  dislocation 

density  between  crossed  polarizers 
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Dislocation  etch  pit  on  (111) 
face  of  Gd3Ga5012  (1000X) 
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would  be  expected.  The  growth  striations  in  the  substrate 
are  reproduced  in  epitaxial  films  grown  by  either  LPE  or 
CVD.  Evidently,  the  surface  of  the  substrate  is  slightly 
attacked  before  film  deposition  begins.  In  some  cases, 
the  striations  in  the  film  are  so  pronounced  that  they 
impede  magnetic  domain  motion.  However,  film  growth  conditions 
can  be  adjusted  to  minimize  the  replication  of  the  striations 
to  the  point  where  they  appear  to  have  no  effect  on 
domain  motion. 

X-ray  diffraction  topography  has  proven  to  be  a  valuable 
tool  to  study  the  growth  striations  as  well  as  other  imper¬ 
fections  in  the  rare-earth  garnets.  The  transmission 
method  developed  by  Lang  (Ref.  25)  is  sketched  in  Fig.  7 
and  is  illustrative  of  the  general  principals  of  x-ray 
topography.  A  narrow  collimated  beam  of  x-rays  is  directed 
upon  the  surface  of  the  crystal.  The  crystal  is  oriented 
so  that  the  beam  is  diffracted  from  a  set  of  planes  with 
Miller  indices  (h  k  1  ) .  In  the  symmetrical  case  the 
diffracting  planes  are  perpendicular  to  the  specimen  faces 
and  the  incident  and  diffracted  beams  make  equal  angles 
with  the  specimen.  A  photographic  plate  is  usually  used  to 
record  the  diffracted  beam  and  a  slit  prevents  the  direct 
beam  from  reaching  the  plate.  The  crystal  and  plate  are 
mounted  so  that  they  can  be  moved  during  exposure  to 
record  a  large  area  of  the  sample.  Bonse,  Hart  and  Newkirk 
(Ref.  26)  have  discussed  other  modes  of  x-ray  topography. 
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Transmission  x-ray  topography  requires  thin  sections 
of  garnet  (about  100pm  thick) .  Usually  it  is  impractical 
to  use  such  thin  samples  for  substrates  for  subsequent 
epitaxial  growth.  Other  variations  of  x-ray  topography 
for  garnets  have  been  evaluated  (Ref.  27)  and  Glass  (Ref.  28) 
has  shown  that  reflection  topography  using  a  double  crystal 
diffractometer  is  especially  useful  for  the  study  of  both 
substrates  and  epitaxial  films.  Reflection  topography 
is  limited,  of  course,  by  the  effective  penetration  of 
x-rays  to  approximately  the  first  20  microns  near  the 
surface  of  the  specimen. 

Much  of  the  information  revealed  by  reflection  topography 
can  also  be  obtained  by  phase  interference  microscopy 
(Ref.  19).  The  phase  interference  microscope  accentuates 
small  surface  irregularities  and  makes  obvious  features 
which  cannot  be  detected  by  ordinary  microscopy.  For 
these  studies,  the  surfaces  should  be  etched  as  described 
the  determination  of  dislocations  above. 

4 • 5  Routine  Substrate  Evaluation 

The  crystal  grower  and  the  epitaxial  growth  technician 
need  a  means  of  rapidly  evaluating  candidate  substrate 
materials.  The  techniques  discussed  here  yield  the  pertinent 
information  required  for  judicious  substrate  selection; 
however,  many  of  the  techniques  are  complementary  and 
some  (x-ray  topography,  for  example)  are  not  really  suitable 
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for  routine  evaluation.  It  is  felt  that  an  adequate  scheme 
for  the  routine  evaluation  of  substrate  materials  would 
involve  the  following. 

Each  boule  should  be  examined  between  crossed-polarizers . 
This  reveals  strain  and  its  general  distribution.  In  some 
cases,  this  examination  alone  can  serve  as  the  basis  for 
the  rejection  of  an  entire  boule. 

Representative  wafers  should  then  be  sliced  from  the 
crystal,  polished  and  etched.  Microscopic  examination  with 
both  transmitted  and  reflected  light  would  reveal  inclusions 
and  dislocations.  These  are  the  most  serious  defects 
commonly  found  in  garnets  and  largely  determine  the  defect 
density  of  the  epitaxial  film. 

Final  examination  of  the  wafers  under  the  phase- 
contrast  microscope  would  reveal  the  extent  of  growth 
striations.  With  this  knowledge,  the  film  grower  can  adjust 
growth  conditions  to  minimize  banding  in  the  epitaxial  film. 
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5.  FILM  THICKNESS 

The  diameter  of  magnetic  bubble  domains  and  the  bias 
field  for  stable  operation  in  a  given  film  are  both  functions 
of  the  film  thickness  (Ref . 29 , 30 , and  31).  Thus,  accurate 
knowledge  of  the  thickness  is  an  important  device  parameter. 
It  is  crucial  in  assembling  a  bubble  domain  system  with 
more  than  a  single  film  that  these  be  matched  in  thickness 
to  within  approximately  1%  as  well  as  uniform  to  that 
tolerance  (Ref. 32) .  in  addition,  many  of  the  magneto- optic 
techniques  to  be  discussed  later  in  this  report  for  the 
determination  of  other  film  properties  require  a  knowledge 
of  film  thickness.  Thus  the  measurement  of  this  elementary 
property  takes  on  a  considerable  importance. 

During  the  stage  in  which  orthoferrites  were  the  most 
interesting  bubble  materials  available  samples  were  generally 
sliced  out  of  flux  grown  bulk  samples.  At  this  stage,  with 
the  entire  thickness  an  active  part  of  the  bubble  domain 
sample,  any  standard  means  for  the  measurement  of  thin 
samples  (such  as  a  micrometer)  was  satisfactory  for  thickness 
determination.  With  the  advent  of  uniaxial  garnet  materials, 
with  their  smaller  bubbles  and  correspondingly  smaller 
thicknesses,  growth  of  films  on  non-magnetic  substrates  was 
investigated.  At  present,  the  epitaxial  growth  of  a 
magnetic  garnet  film  on  non-magnetic  garnet  substrate  is 
by  far  the  most  widely  used  means  for  preparing  magnetic 
bubble  samples. 


5.1  Edge  Viewing  Methods 


The  first  method  for  the  determination  of  film  thickness 
in  these  samples  probably  involved  breaking  the  sample  and 
observing  the  cracked  edge  microscopically.  The  magnetic 
film  can  be  distinguished  from  the  substrate  either  by  its 
greater  optical  absorption  (in  transmitted  light)  or  by 
etching  in  hot  phosphoric  acid  to  reveal  the  boundary.  This 
method  is  useful  on  rough  deposits  and  requires  less 
equipment  than  those  outlined  below  but  has  the  obvious 
disadvantages  of  being  destructive  to  the  sample  and  of 
limited  accuracy.  In  addition  the  garnets  do  not  cleave 
well  so  that  one  must  usually  deal  with  a  jagged  edge  and 
hope  for  a  sufficiently  flat  region  to  reveal  the  boundary 
sought.  Such  measurements,  because  of  their  directness, 
have  been  useful  in  checking  some  of  the  other  techniques 
to  be  discussed. 

A  more  reproducible  and  accurate  varient  of  the  simple 
break  approach  is  metallurgical  cross-sectioning  (Ref.  33  and 
34) .  Because  of  the  hardness  of  the  garnets  it  is  difficult 
to  find  a  material  which  is  sufficiently  hard  to  avoid 
faster  removal  of  the  material  in  which  the  sample  is 
potted  and  the  rounding  of  the  sample's  edges  which  follows. 
This  effect  poses  a  limitation  of  several  percent  on  the 
accuracy  with  which  the  thickness  of  a  typical  film  can  be 
measured  in  this  way.  In  addition  the  method  remains 
essentially  destructive. 
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5. 2  Interferometric  Methods 

As  film  uniformity  improved  non-destructive  techniques 

based  upon  optical  interferometry  became  possible.  In  these, 

a  beam  of  monochromatic  light  is  directed  onto  the  sample 

and  either  the  transmitted  or  reflected  beam  intensity 

measured  as  the  wavelength  of  the  incident  light  is  varied. 

Required  for  the  technique  are  a  flat  sample  surface,  good 

film  uniformity  over  the  region  sampled,  and  a  significant 

and  sharply  defined  difference  in  the  index  of  refraction 

between  the  film  and  the  substrate.  The  fraction  of  the 

incident  light  reflected  at  the  film  (index  nf)  -  substrate 

(index  n  )  interface  is  given  by 
s 


Since  this  reflected  light  forms  an  essential  part  of  the 
interference  phenomenon  the  requirement  for  finite  nf-n 
arises.  Fortunately  the  iron  ions  which  are  present  in  the 
film  but  not  in  the  substrates  of  typical  garnet  samples 
yield  considerable  increase  in  index  as  well  as  absorption 
(see  Section  17).  Typically,  n  ~  1.70  (Pef.  35)  while 

5 

nf  ~  2.35  in  the  red  portion  of  the  spectrum  (Ref.  36)  so 
that  R2  =  0.026  which  is  quite  adequate  to  observe  inter¬ 
ference  effects. 

The  basic  interference  effect  for  thin  films  is 
treated  in  texts  on  optics  (Ref .37  and  38)  and  will  not 
be  repeated  here.  The  optical  path  difference  between  the 
light  reflected  from  the  free  surface  of  the  film  and  the 
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interface  with  the  substrate  causes  various  phase  differences 
to  occur  as  X,  the  incident  light  wavelength,  is  varied. 
These,  in  turn,  cause  the  oscillations  in  intensity  typical 
of  the  phenomenon.  The  effect  has  been  employed  in  the 
garnets  both  in  transmission  (Ref.  35)  and  reflection 
(Refs.  36,  39, and  40).  Reflection  is  to  be  preferred  for 
the  following  reasons:  1)  The  back  surface  of  the 
substrate  is  unimportant  in  reflection  while  an  optical 
polish  is  desirable  for  the  transmission  measurement. 

Thus,  the  reflection  technique  requires  less  sample 
preparation.  2)  The  oscillations  are  a  larger  fraction 
of  the  total  signal  in  reflection.  This  can  be  seen  as 
follows:  The  fraction  of  the  total  light  reflected  at  the 
front  film  surface  is 


(nf  -  l)2 
(nf  +  1)* 


0.16 


(5-2) 


Neglecting  secondary  internal  reflections  and  absorption 
in  the  film  this  fraction  interferes  with  the  0.026  fraction 
reflected  at  the  film- substrate  interface.  Taking  into 
account  the  diminution  of  this  internal  beam  by  the  first 
surface  reflection,  the  internally  reflected  beam  is  ~  14% 
of  the  front  surface  reflected  beam.  This  leads  to  a 
maximum  fractional  variation  of  the  reflected  beam  due  to 
interference  of  ~  28%.  Since  the  transmitted  beam,  after 
consideration  of  reflection  at  the  backside  of  the  substrate, 
is  approximately  70%  of  the  incident  beam,  the  maximum 


34 


fractional  variation  of  this  beam  will  be  ~  6.4%.  Those 
considerations  are  qualitatively  borne  out  by  comparison 
of  the  results  (Ref.  35)  in  transmission  and  (Ref.  39) 
in  reflection.  (See  also  Fic.  8a) 

The  calculation  of  the  film  thickness  from  these  data 
follows  from  a  consideration  of  the  extra  path  length  in  the 
film.  For  normal  incidence  on  a  film  of  index  of  refraction 
n,  this  yields 


h  = 


(5-3) 


where  A^+^  and  A^  are  wavelengths  corresponding  to 
successive  maxima  or  minima.  In  practice  an  average 
spacing  in  1/A  is  calculated  as 


(5-4) 


in  order  to  improve  accuracy.  The  use  of  this  equation 
assumes  a  constant  (or  properly  weighted  average)  index, 
n  (an  assumption  which  can  be  avoided) . 

Roman,  et  al  (Ref.  36)  have  assumed  a  linear  variation 
of  n  with  X  so  that  the  index,  n(X2)  at  X^  is  given  in  terms 
of  n(Xj)  at  X^  and  the  variation  of  n  with  X  as 

n(X2)  =  nfA^  +  ^X2  -  A,)  (5-5) 


Then  h 


N  f  Vl  A 

2  \X2~XlJ 


n  (X1) 


(5-6) 
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where  N  is  the  number  of  intensity  oscillations  between 

X  and  X_.  As  those  authors  point  out  this  equation  can 
1  2 

be  used  as  a  guide  to  choosing  the  proper  weighted  average 
for  use  in  the  simpler  equation  5_3 .  Within  the  visible 
to  near-infrared  region  gy  <  0  (Ref .41  and  42)  so  that  an 
index  somewhat  larger  than  those  corresponding  to  the  wave¬ 
length  range  used  should  be  employed.  While  this  conclusion 
appears  surprising,  because  of  the  nature  of  the  equation 
it  is,  in  fact,  correct  as  can  be  seen  by  direct  substitution 
in  equation  5—7.  Roman  et  al  quote  values  of 
index  for  three  bubble  materials  as  follows t 

O 

Material  n  at  X  =  6330A 


Gd2.3Tb0.7Fe5°12 
EulEr2Ga0 . 7Fe4 . 3° 12 
Er2GdlGa0.7Fe4.3°12 


2.42  +  0.05 

2.35  +  0.05 

2.36  +  0.05 


The  value  of  Sv  near  6000A  can  be  inferred  from  unpublished 
dA 

values  of  Tabor  as  quoted  by  Shumate  (Ref.  40) .  These  are 


X  n 

6500A  2.350 

5 500 A  2.396 

leading  to  =  -4.6  x  10‘5/A.  Using  this  a  reasonable 
estimate  of  n  for  equation  5-3  and  the  convenient  wavelength 
range  6,000  to  10,000A  is  2.45. 

Shumate  (Ref.  40)  has  proposed  the  equation 
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_  N  2.396  2.350  ~1 

2  X1  X2 


(5-7) 


where  *2  and  Xj  are  wavelengths  corresponding  to  extrema 
of  the  same  sign  near  5500A  and  6500A  respectively  and  N  is 
the  number  of  intensity  oscillations  between  these  wavelengths. 
This  is  the  most  accurate  of  the  equations  quoted  here  in  that 
the  indices  used  are  those  measured  at  or  near  the  wave¬ 
lengths  employed.  It  also  requires  the  best  films  since  small 
non-uniformities  reduce  the  oscillations  at  short  wave¬ 
lengths  most  strongly.  Such  an  approach,  perhaps  extended 
to  somewhat  longer  wavelengths  when  the  full  dispersion 
curve  is  known,  will  probably  be  the  most  accurate  and 
convenient  means  of  determining  bubble  film  thickness 
for  the  foreseeable  future. 


Systems  for  the  measurement  of  such  interference 
traces  are  available  from  many  suppliers  of  monochromators. 
A  relatively  inexpensive  committed  system,  together  with 
typical  data,  is  shown  in  Figure  8  a  &  b  as  assembled  at 
Monsanto.  Specific  apparatus  used  has  been  listed  pre¬ 
viously  (Ref.  39).  when  used,  over  the  full  range,  the 
system  can  approximate  the  relative  accuracy  of  1% 
necessary  to  match  films  in  multi-film  devices.  Absolute 
accuracy  is  +  4%,  however,  as  a  result  of  uncertainty 
in  n  and  beam  divergence  (normal  incidence  is  assumed  in  all 
of  the  above  equations).  Variations  in  thickness  from  one 
point  to  another  on  a  sample  can  be  handled  as  described 
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Figure  8b 

Figure  8.  Thickness  system  and  typical  data  curve 
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in  the  next  section. 

The  system  in  use  at  Bell  Telephone  Laboratories  (Ref.  40) 
differs  from  this  in  that  the  sample  is  visually  observed 
(by  eye  in  the  visible  or  with  the  aid  of  a  silicon  vidicon 
tube  in  the  infrared)  and  the  number  of  interference 
fringes  passing  a  given  point  as  the  wavelength  is  varied 
in  simply  counted.  This  system  has  the  advantage  of  pro¬ 
viding  both  thickness  and  thickness  variation  information 
simultaneously  but  requires  the  investment  in  a  closed 
circuit  TV  system  for  the  infrared  region.  In  addition, 
as  film  uniformity  approaches  the  desired  flatness, 
band  fringes  are  no  longer  visible  and  the  entire  surface 
oscillates  in  reflectivity  as  the  wavelength  changes. 

In  this  case  the  quantitative  data  provided  by  a  normalized 
photomultiplier  or  solar  cell  system  is  a  necessity  in 
order  to  find  the  extrema  accurately.  It  is  suggested 
that  a  system  utilizing  photoelectric  readout  is  to  be 
preferred  and  the  separate  but  simple  system  of  the  next 
section  used  to  determine  thickness  variations. 


6.  THICKNESS  VARIATIONS 

The  observation  of  a  broken  or  cut  and  polished  edge 
can  yield  information  on  thickness  variations  along  that 
boundary  but  these  destructive  techniques  are  undesirable 
as  production  procedures.  Fortunately,  optical  interference 
again  provides  a  very  convenient  alternative.  Interference 
effects  entirely  similar  to  those  described  in  Section  5 
when  wavelength  is  varied  would  also  be  observed  if  the 
film  thickness  were  varied  at  constant  wavelength. 

Such  an  effect  is  actually  observed  when  regions  of 
different  thickness  are  present  in  the  same  sample.  In 
this  case  the  dark  bands  arising  from  destructive 
interference  in  the  reflected  light  correspond  to  contours 
of  constant  film  thickness  satisfying  the  relationship 
2  h  n  =  ml.  Once  again  h  is  the  film  thickness,  n,  the 
index  of  refraction  at  wavelength  A,  and  m,  an  integer. 

The  difference  in  thickness  between  successive  dark 
bands  is  then  Ah  =  A/2n. 

The  workers  at  Bell  Telephone  Laboratories  have 
generally  observed  the  fringes  in  the  infrared  region  using 
monochromator  illumination  and  a  silicon  vidicon  tube 
(Ref.  36).  This  system  permits  simultaneous  determination 
of  thickness  (as  mentioned  in  the  previous  section)  and  its 
variation.  While  some  photographs  of  these  fringe  patterns 
have  been  made  the  usual  system  used  at  BTL  for  generating 
a  map  of  thickness  variations  has  been  to  mount  the  sample 
on  a  microscope  stage  equipped  with  position  transducers, 


illuminate  it  with  monochromatic  light,  and  generate  a 

map  on  a  X-Y  recorder  by  manually  tracing  the  fringes 
(Ref.  40). 

A  more  rapid  and  convenient  system  is  in  use  at 
Monsanto  (Ref.  39).  This  consists  simply  of  taking  a 
close-up  photograph  of  the  sample  illuminated  by  mono¬ 
chromatic  light  next  to  a  grid  using  the  system  of 
Figure  9a.  The  helium  filled  discharge  lamp  provides  a 
strong  line  at  0.5876ym  wavelength.  This  yields  a 
thickness  difference  between  the  successive  dark  fringes 
of  (0.122  +  0.005) urn  based  upon  n  =  2.40  at  this  wavelength. 
These  fringes  are  readily  visible  to  the  eye  when  viewed 
in  this  light.  This  technique  has  proved  invaluable 
in  our  efforts  to  improve  film  flatness.  In  order  to 
photograph  the  interference  fringes  a  glass  beam  splitter 
is  used  as  indicated  in  the  figure.  A  typical  result 
is  shown  in  Figure  9b. 

Such  a  map  in  itself  does  not  show  the  sign  of  the 
thickness  gradient.  This  can  often  be  determined  by  visual 
inspection  since  a  significantly  thicker  region  will 
appear  darker  green.  If  the  thickness  variations  are  too 
small  to  yield  detectable  color  differences  another  simple 
technique  is  still  available.  This  involves  rotating  the 
sample  while  viewing  the  fringes  under  monochromatic  light. 

As  the  sample  is  rotated  away  from  a  perpendicular  viewing 
angle  the  dark  bands  appear  to  move  toward  thicker  parts 
of  the  film.  This  can  be  seen  as  follows:  The  optical 


Figure  9b 


Figure  9.  Thickness  variation  system  and 
typical  thickness  map 


42 


path  difference  between  the  light  reflected  from  the 
front  surface  and  that  reflected  from  the  film-substrate 
interface  is  given  by  A  (OP)  =  2  n  h  cos  0'  where  0'  is 
the  angle  between  the  normal  to  the  film  and  the  propagation 
direction  of  the  light  internal  to  the  film  (Ref.  37).  As 
the  sample  is  rotated  away  from  the  perpendicular,  0' 
increases  from  zero  and  the  optical  path  difference  at  the 
spot  of  thickness  h  decreases.  However,  the  original 
dark  band  will  move  to  a  point  at  which  A  (OP)  is  the  same 
as  it  originally  was  at  the  spot  of  thickness  h,  i.e.  it 
must  move  in  the  direction  of  increasing  h. 

The  combination  of  variable  wavelength  interferometry, 
monochromatic  light  photography,  and  visual  inspection  as 
a  function  of  angle  under  the  monochromatic  lighi: 
therefore  provides  complete  thickness  information  at  all 
points  of  the  films.  Other  monochromatic  light  sources, 
such  as  a  sodium  lamp  or  suitably  diffused  laser  beam, 
are  satisfactory  for  viewing  the  fringes.  The  shortest 
possible  wavelength  is  to  be  preferred,  however,  because 
of  its  greater  sensitivity  to  small  thickness  variations. 

O 

For  example,  using  the  5500A  index  quoted  by  Shumate, 
successive  fringes  at  that  wavelength  correspond  to  0.115pm 
thickness  difference,  to  be  compared  with  the  0.122pm  figure 

O 

quoted  above  for  ~  5900A.  Because  of  the  rapidly  increasing 

O 

iron  absorption,  5000A  is  an  effective  lower  wavelength 
limit  for  the  technique. 
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7.  DEFECT  DETECTION  AND  LOCATION 

Defects  in  a  magnetic  bubble  material  which  impede  the 
motion  of  bubble  domains  will,  at  the  very  least,  reduce  device 
operating  margins  and,  the  worst  case,  may  render  the  sample 
useless  for  device  operation.  Detection  and  location  of  such 
defects  thus  becomes  an  important  phase  of  quality  control, 
both  as  a  check  on  damage  introduced  at  various  stages  of 
sample  processing  and  for  the  selection  of  device  quality  areas. 

Several  of  the  standard  methods  for  the  study  of  defects 
are  applicable  to  magnetic  bubble  materials.  These  include 
etching  studies  (Reference  24  and  43  )  X-ray  topography 
(Reference  44  and  45  )  and  microscopic  examination  (Reference  46  ) . 
These  techniques  are  essentially  the  same  as  those  used  in  the 
evaluation  of  substrates  and  have  been  discussed  in  Section  4. 

We  confine  our  attention  below  to  those  non-destructive 
techniques  uniquely  suited  to  magnetic  bubble  materials.  All 
of  these  involve  the  microscopic  examination  of  films  or  thin 
slabs  in  transmitted,  polarized  light.  The  plane  of 
polarization  is  rotated  by  the  magnetic  material  in  a  sense 
dependent  upon  its  magnetization  direction  (the  Faraday  effect, 
see  Section  17)  and  the  domain  pattern  is  revealed  by  placing 
another  polarizer  (the  analyzer)  in  the  emerging  light  beam  so 
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as  to  transmit  nore  of  the  light  passing  through  one  sign  of 
domain  than  the  other.  The  interaction  of  domain  walls  with 
crystalline  defects  can  then  be  revealed  in  a  number  of  possible 
configurations . 

In  the  orthoferrites  it  is  fairly  straightforward  to 

stabilize  a  single  straight  domain  wall  by  locating  the  sample 

in  a  region  in  which  the  field  strength  perpendicular  to,  the 

sample  plane  is  varying  sufficiently  rapidly  with  position,  i.e., 

dH/dx  >  0  with  H  directed  along  the  z  axis  normal  to  the 
c 

sample  plane.  Hagedorn  (Reference  47  )  has  calculated  0c  in 

terms  of  the  material  saturation  magnetization,  4ttm  ,  and 

s 

2  2 

characteristic  length, £  (5  aw/4TTMs  where  ow=wall  energy/cm  )  ,  as 

4TTM 

0  =  0.016  .  7-1 

C  Jv 

A  very  versatile  system  for  the  generation  of  this  type  of 
field  configuration  using  two  permanent  magnets  has  been 
described  by  Shumate  (Reference  47)  and  is  reproduced  here 
(Figure  10).  When  the  sample  is  now  moved  within  this  gradient 
field  and  the  domain  wall  observed  microscopically,  defects 
which  interact  with  the  wall  will  cause  distortions  from  the 
straight  configuration,  i.e.,  from  the  line  corresponding  to 
zero  applied  field.  Examples  of  such  distortions  have  been 
shown  by  Burmeister,  et  al  (Reference  48  and  49).  The  same 
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Figure  10.  Device  for  establishing  variable  gradient  field. 

After  Shumate,  ref.  47.  The  sample  rests  on  a 
slide  in  the  gap  between  the  magnets. 


technique  has  been  used  by  Kurtzig  to  investigate  the  inter¬ 
action  of  a  domain  wall  with  specific  defects  such  as  twin 
and  grain  boundaries  (Reference  50)  and  dislocations  (Reference 
44).  If  the  magnitude  of  the  field  gradient  is  known  it  is 
possible  to  calculate  the  coercivity,  Hc»  associated  with  the 

individual  defect  as  H  =  Ax(dH/cb:)  where  Ax  is  the  displacement 

c 

of  the  domain  wall  from  its  equilibrium  position  in  the  gradient 
field  at  the  time  the  vail  breaks  free  of  the  defect.  For  the 
cases  investigated  by  Kurtzig,  this  coercivity  was  found  to  be 
highly  anisotropic,  as  might  be  expected.  Further  discussion 
of  defects  detected  in  this  way  has  been  presented  by  Nielsen 
(Reference  51) . 

The  generation  of  a  single  straight  domain  wall  is 

considerably  more  difficult  in  the  garnets  than  it  is  in  the 

orthoferrites.  The  garnets  are  of  practical  interest  because 

of  their  smaller  bubble  size  relative  to  the  orthoferrites. 

This  bubble  size,  in  turn,  results  from  lower  domain  wall 

energy  and  higher  magnetization.  For  example,  a  typical  garnet 

2 

has  a  ~  0.2  ergs/cm  and  4ttM  ~  150G  while  typical  orthoferrite 
w  3  s 

2 

values  are  o  ~  1  erg/cm  and  4ttM  ~  100G.  The  resulting 
w  5  s 

increase  in  5  required  to  generate  the  single  straight  wall 
c 

is  nearly  a  factor  of  20  and  puts  the  straight  wall  technique 
out  of  practical  reach  for  many  of  the  garnets.  In  addition. 


as  the  required  gradient  becomes  larger,  the  sensitivity 
(Ax  caused  by  a  defect  of  given  pinning  strength)  goes  down 
and  with  it  the  usefulness  of  the  method. 

In  the  presence  of  a  gradient  insufficient  to  bring  about 
a  straight  wall,  a  single  wall  is  still  possible.  This  wall 
will  take  on  the  wavy  to  finger-like  forms  shown  by  Hagedorn 
(Reference  52)  depending  on  the  strength  of  the  gradient 
applied.  The  passage  of  such  a  pattern  across  a  magnetic 
bubble  film  is  still  sensitive  to  the  presence  of  defects. 

This  technique  is  being  used  for  the  detection  of  defects  at 
Bell  Northern  of  Canada  where  it  is  thought  to  be  more 
sensitive  than  the  methods  to  be  described  below  (Reference  53)- 

Another  method  for  defect  detection  particularly  suitable 
for  the  orthoferrites  has  been  developed  by  Fischer  and 
Shumate  (Reference  54).  In  it  a  bubble  domain  is  moved  across 
the  sample  by  means  of  a  magnetized  probe  and  distortions  of 
the  bubble  caused  by  defects  are  detected  optically.  This 
requires  a  bias  field  to  be  applied  opposite  to  the  magnetization 
in  the  bubble  and  the  probe  which  sometimes  causes  reversal  of 
the  probe  magnetization  and  a  temporary  stop  to  the  scan 
(Reference  47).  The  raster  scanning  mechanism  which  was  built 
to  automate  this  system  makes  it  a  possible  production  technique 


48 


for  the  orthoferrites.  As  bubble  size  becomes  smaller, 
however,  the  number  of  passes  required  in  the  raster  increases 
and  the  photoelectric  signal  involved  in  a  bubble  distortion 
decreases  making  the  technique  marginal  in  the  garnets. 

Nikitenko,  et  al  (Reference  55)  have  reported  studies  of 
the  influence  of  dislocations  on  magnetic  domains  in  gallium 
substituted  YIG.  The  domains  were  observed  by  the  Faraday 
effect  in  the  usual  fashion  while  the  dislocations  were 
simultaneously  observable  due  to  the  birefringence  brought 
about  by  their  strain  fields.  They  were  able  to  observe  the 
nucleation  and  binding  of  domains  to  the  dislocations  but  did 
not  determine  quantitative  binding  energies.  Their  samples 
were  quite  thick  (0.3  to  2mm)  and  observations  made  at  1.2|-t 
wavelength.  It  seems  unlikely  that  as  clear  cut  defect 
observations  can  be  made  by  this  technique  using  the  much 
thinner  samples  of  interest  in  practical  bubble  devices. 

The  most  widely  used  techniques  for  defect  detection  in 
the  garnets  are  based  upon  the  observation  of  domain  motion 
in  an  alternating  bias  field.  As  described  by  Geusic,  et  al 

(Reference  56)  the  bias  field  is  cycled  slowly  with  time  and 
points  which  interfere  with  the  motion  of  strip  domains  are 
detected  visually.  Two  photomicrograph  examples  are  presented 
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in  that  paper.  Shumate,  in  reporting  on  this  method 
(Reference  47) ,  indicated  that  a  60  Hz  sine  wave  bias  is 
adequate,  although  judicious  addition  of  a  dc  bias  as  well 
often  aids  in  elucidating  weak  imperfections  and  gradients. 

In  this  60  Hz  technique  a  sufficient  AC  field  amplitude  is 
applied  to  cause  appreciable  domain  motion  when  viewed  in  a 
polarizing  microscope.  Regions  in  which  domains  are  still 
pinned  at  defects  will  be  decorated  by  non-moving  walls.  The 
technique  requires  a  judicious  choice  of  field  amplitude, 
since  a  large  amplitude  will  tear  domains  away  from  weak 
pinning  sites  and  those  sites  will  go  undetected.  In  practice 
the  minimum  amplitude  which  causes  a  blurred  condition  serves 
very  well  to  delineate  strong  defects  or  those  with  some 

spatial  extent  such  as  scratches  and  clusters  of  point  defects. 
For  weaker  point  defects  it  is  often  necessary  to  go  to  lower 
amplitudes,  add  DC  bias,  and  look  for  points  at  which  the 
domain  pattern  tends  not  to  move  or  to  move  in  a  "snap-action" 
fashion.  Hiskes  (Reference  57)  and  Shumate  (Reference  40)  have 
presented  photomicrographic  examples  of  this  technique.  We 
show  an  additional  example  of  a  line  defect  in  Figure  11. 

Chaudhari  (Reference  18)  has  recently  described  another 
magneto-optic  technique  for  the  detection  of  defects.  An  array 
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Figure  11a.  Domain 


pattern  without  bias  field 


Figure 

Figure 


ing  field 

11.  Defect  detection  in  oscillating  bias  field 


of  bubbles  is  made  to  flow  as  a  raft  (and  presumably  generated 
as  well)  by  a  suitably  formed  pulse  coil.  Defects  are  detected 
by  hang  up  of  individual  bubbles  or  by  the  distortion  of  the 
raft  from  hexagonal  symmetry.  Several  examples  and  the 
arguments  in  favor  of  this  approach  have  been  presented  by 
Chaudhari.  Unfortunately  no  detailed  comparisons  as  to 
sensitivity  or  visibility  were  given.  We  feel  that,  while  it 
may  be  too  early  for  a  final  decision,  the  oscillating  stripe 
or  moving  convoluted  wall  both  provide  more  useful  techniques 
for  the  routine  scanning  of  samples  for  defects. 

Defect  maps  are  readily  made  by  connecting  the  microscope 
stage  to  position  transducers  which,  in  turn,  communicate 
with  an  X-Y  recorder.  The  system  used  at  Monsanto  and  a 
typical  defect  map  are  shown  in  Reference  39.  The  sample  is 
manually  scanned  at  200X  by  an  operator  who,  upon  observing 
a  region  of  domain  wall  pinning,  centers  it  in  the  microscope 
field  of  view  and  presses  a  footswitch  to  drop  the  recorder 
pen  and  mark  the  map.  Sixty  (60)  cycle  modulation  derived  from 
a  variable  transformer  is  quite  satisfactory. 

It  is  difficult  at  this  point  to  say,  in  any  of  these 
techniques,  whether  all  significant  defects  are  being  detected 
or  whether  some  of  those  detected  may,  in  fact,  be  too  weak  to 
cause  problems.  What  is  needed  in  this  regard  is  a  correlated 
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study  of  defect  strength  from  the  above  measurements  and 
bubble  propagation  errors  under  a  standard  circuit  mask.  The 
publication  of  a  study  of  defects  using  both  the  convoluted 
single  boundary  stabilized  by  a  gradient  and  the  modulated 
strip  domain  technique  on  the  same  specimens  would  also  be  a 
service.  At  the  present  time  we  recommend  the  use  of  either 
of  these  approaches  for  the  garnets. 


8. 


CHARACTERISTIC  LENGTH  AND  DOMAIN  DIMENSIONS 


The  characteristic  length,  l,  to  be  used  in  the  present 

report  follows  the  notation  introduced  by  Thiele  (Reference  29 

and  30)  and  now  rather  widely  used  in  the  literature  of 

2 

cylindrical  domain  technology,  namely,  SL  =  o  /4ttm  .  Here  a 

w  s  w 

is  the  energy  per  unit  area  of  a  domain  wall  and  M  is  the 

s 

saturation  magnetization.  This  parameter  is  clearly  a 
combination  of  two  more  basic  ones.  It  is  of  interest, 
hovever,  because  it  governs  the  scale  of  domain  sizes 
characteristic  of  a  given  magnetic  bubble  material.  Conversely 
the  measurement  of  certain  domain  sizes  and  the  sample  thickness 
is  sufficient  to  yield  the  value  for  1  so  that  only  one  further 
determination  (generally  4ttm  )  is  required  to  specify  both  a 

S  w 

and  4ttMs.  In  this  report  we  will  devote  a  separate  section  to 
because  there  exist  a  number  of  alternative  approaches  to 
its  measurement  which  are  worthy  of  note. 

8.1  Methods  Employing  Bubbles 

The  theory  of  Thiele  begins  with  the  assumption  of  three 
independent  contributions  to  the  total  sample  energy  for  a 
sample  containing  a  single  domain,  the  energy  due  to  domain 
walls,  that  due  to  interaction  of  the  sample  with  the  external 
field,  and  that  due  to  internal  magnetostatic  interactions. 
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The  solutions  for  the  limits  of  stability  of  a  cylindrical 
magnetic  domain  are  presented  in  graphical  form  and  reproduced 
here  as  Figure  12.  Here  SQ  (d/h)  represents  the  stability 
limit  against  collapse  and  S2  (d/h)  the  stability  limit  against 
distortions  of  twofold  symmetry  (i.e.,  strip-out).  Thus, 
cylindrical  domains  are  stable  in  the  area  bounded  by  Sq  and 
S^.  The  ordinate  for  these  curves  corresponds  to  the  ratio 
H/h  (where  h  is  the  sample  thickness)  and 
therefore  forms  the  basis  for  the  determination  of  £  from 
measurements  of  bubble  diameters  and  sample  thickness.  For 
this  simplest  use  of  the  Thiele  theory  the  diameter  either  at 
the  strip  to  bubble  transition  or  at  collapse  is  required. 

Several  difficulties  are  encountered  in  carrying  out  this 
approach.  1)  The  accuracy  of  the  strip  to  bubble  transition 
point  is  limited  both  by  coercivity  effects  in  real  samples 
and  by  the  ability  of  the  eye  to  detect  small  non-circular 
distortions.  2)  Once  collapse  of  the  bubble  of  interest  is 
established  it  is  no  longer  available  for  measurement.  Thus 
one  must  either  shift  to  a  nearby  (but  still  "isolated”)  bubble 
for  the  diameter  determination  or  take  a  series  of  magnetic 
field-diameter  data  as  collapse  is  approached  and  extrapolate 
the  resulting  curve  to  the  collapse  diameter  using  the  measured 
field  at  collapse.  3)  The  measurement  of  domain  diameters. 


Magnetostatic  Force-Stability 


Domain  Diameter  to  Thickness  Ratio  -  d/h 


Figure  12.  The  stability  curves  for  collapse  (SQ)  and 
strip  out  (S2)  for  bubble  domains  and  the 
radial  force  function  F. 

After  Thiele,  ref.  30. 
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while  relatively  simple  and  accurate  in  the  ortho^errites 
with  their  large  domains  (>  25um) ,  becomes  less  accurate  for 
the  smaller  domains  of  the  garnets.  Hagedorn,  et  al 
(Reference  58) ,  in  a  test  of  the  accuracy  which  the  Thiele 
approach  permits,  quote  an  uncertainty  in  bubble  diameter  of 
+,  0.2|im  which,  when  coupled  to  the  curve  fitting  approach  to 
be  outlined  next,  leads  to  a  relative  accuracy  in  l  of  +  7%. 

The  experience  at  Monsanto  indicates  that  a  diameter 
uncertainty  of  +  0.5um  is  realistic  for  good  but  not  ultimate 
optical  equipment.  The  resulting  uncertainty  in  is  not 
satisfactory  in  light  of  the  alternatives  to  be  discussed 
below. 

The  curve  fitting  approach  of  Hagedorn,  et  al  (Reference  58) 
involves  making  a  least-squares  fit  of  a  series  of  field  and 
corresponding  diameter  values  for  a  film  of  known  thickness 
to  the  equilibrium  equation  derived  by  Thiele  (Reference  29, 
Equation  69) . 
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relative  +  7%  accuracy  quoted  above.  This  mechod  does  not 

hold  great  promise  as  a  production  line  technique  of  reasonable 
accuracy. 

8.2  Isolated  Straight  Wall  Method 

A  rather  unique  approach  to  the  determination  of  a  (and, 
with  less  accuracy,  4ttms)  has  been  outlined  by  Hagedorn 
(Reference  52).  It  makes  use  of  the  isolated  straight  wall 
configuration  mentioned  in  the  previous  section.  As  the  field 
gradient  needed  to  stabilize  this  single  wall  is  reduced, 
deviations  from  straightness  develop  which  have  a  periodicity 
and  critical  field  gradient  characteristic  of  the  sample 
thickness  and  the  material  under  test.  The  characteristic 
length  is  a  function  only  of  the  periodicity  of  the  distortion 
and  the  thickness  near  the  onset  of  this  phenomenon  and  can  be 
determined  with  good  accuracy.  4TTMg  requires  a  knowledge  of 
the  gradient  at  which  the  distortion  sets  in  and  this  is  not 
so  well  defined.  Unfortunately,  as  mentioned  previously,  the 
generation  of  an  isolated,  straight  wall  in  the  garnets  is  a 
difficult  project,  requiring  very  large  gradients,  so  that 
this  approach  has  not  been  reported  for  these  materials. 

8*3  Methods  Employing  Strip  Domains 

The  normal  domain  pattern  for  a  magnetic  bubble  sample 
in  zero  applied  field  is  a  serpentine  array  of  strips  (see. 


for  example.  Figure  11).  The  alternate  regions  of  magnetization 

into  and  out  of  the  sample  will  have  equal  width  in  order  to 

yield  a  net  magnetic  moment  of  zero.  This  width  is  a  function 

only  of  the  sample  thickness  and  the  characteristic  length  in 

the  limit  of  large  anisotropy  assumed  by  the  theories.  Fowlis 

and  Copeland  (Reference  59)  and  Shaw,  et  al  (Reference  60  and 

3l)  have  solved  the  earlier  theories  of  Malek  and  Kambersky 

(Reference  61)  and  Kooy  and  Enz  (Reference  62)  to  yield  fc/h 

from  measured  values  of  P  and  h.  Here  P  is  the  spatial 

period  or  repeat  distance  of  a  set  of  parallel  strip  domains 

measured  in  the  direction  perpendicular  to  the  strips.  The 

results  of  these  calculations  are  shown  in  Figure  13 

and  are  presented  in  tabular  form  in  Table  1. 

As  mentioned  above,  the  equations  upon  which  these  results 

are  based  make  the  assumption  of  high  anisotropy,  i.e.,  H  » 

A 

4rrMg.  Here  Ha  is  the  uniaxial  anisotropy  field  which  is  related 

to  the  frequently  used  anisotropy  constant  K  as  H  =  2Ku  . 

u  A  M 

s 

The  above  inequality  is  not  fully  achieved  in  practical  materials 
nor  is  it,  in  fact,  desirable.  As  Hagedorn  and  others 
(References  63,  64,  and  65)  have  pointed  out  the  domain  wall 
mobility  is  inversely  proportional  to  Ku  so  that  a  high  aniso- 
tropy  implies  low  mobility,  quite  undesirable  in  a  magnetic 


59 


m 


mmmmHi*, 


TABLE  I 

Characteristic  length.  Si,  divided  by  thickness,  h,  vs.  strip 
domain  period,  PQ,  divided  by  thickness. 


Vh 

SL/h 

Vh 

Jl/h 

0.50 

.0085 

6.20 

.389 

0.60 

.0122 

6.40 

.398 

0 . 70 

.0166 

6.60 

.406 

0.80 

.0216 

6.80 

.415 

0.90 

.0273 

7.00 

.424 

1.00 

.0335 

7.20 

.433 

1.20 

.0473 

7.40 

.441 

1.40 

.0626 

7.60 

.448 

1.60 

.0788 

7.80 

.457 

1.80 

.0956 

8.00 

.465 

2.00 

.1126 

8.20 

.472 

2.20 

.1295 

8.40 

.478 

2.40 

.1464 

8.60 

.485 

2.60 

.1629 

8.80 

.493 

2.80 

.1789 

9.00 

.500 

3.00 

.1945 

9.20 

.508 

3.20 

.2100 

9.40 

.  514 

3.40 

.2247 

9.60 

.521 

3.60 

.2387 

9.80 

.526 

3.80 

.2525 

10.00 

.533 

4.00 

.2660 

10.50 

.546 

4.20 

.2790 

11.00 

.562 

4.40 

.2915 

11,50 

.575 

4.60 

.304 

12.00 

.588 

4.80 

.315 

12.50 

.602 

5.00 

.327 

13.00 

.613 

5.20 

.338 

13.50 

.625 

5.40 

.348 

14.00 

.637 

5.60 

.360 

14.50 

.648 

5.80 

.370 

15.00 

.659 

6.00 

.379 
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bubble  device  material.  The  Thiele  theory  for  bubble  domains 
was  also  based  upon  this  assumption  of  large  anisotropy.  In 
this  case  the  recent  results  of  Lin  and  Tu  (Reference  66) 
and  particularly  those  of  DeBonte  (Reference  67)  extend  the 
theory  to  finite  anisotropy.  At  Monsanto  we  have  carried 
out  more  accurate  calculations  including  finite  anisotropy 
effects  for  strip  domains  based  upon  the  Equations  1  through 
3  of  Kooy  and  Enz.  The  equilibrium  equation  which  results  is 


h 


00 


n=l 


sinh(nrra)  fsinh(nm)  +  Jy  cosh(nTTg)  -nna  /u  1 
[sinh(nm)+  1  cosh(nm)  ]  2 


n  odd  (8-2) 

Here  0=1  +  4ttMs/Ha  and  a  =  ^-  .  (Ref.  68) 

o 

The  results  of  this  calculation  for  values  of  p  =  1.0, 
1.25,  and  1.5  (i.e.,  H  /4ttm  =  °°,  4,  and  2)  are  given  in 
Table  2  .  Note  that  the  total  variation  of  S/h  near  the 
practical,  interesting  ranqe  2,/h  ~  0.25  (for  example,  at 
PQ/h  =4.0)  is  only  4%.  Nevertheless,  as  measurements  become 
more  accurate  such  an  effect  should  be  taken  into  account. 

Essentially  all  domain  size  measurements  are  presently 


being  done  in  polarizing  microscopes  equipped  with  a  graduated 
eyepiece.  The  entire  optical  system  is  calibrated  by  placing 
an  accurately  ruled  slide  on  the  microscope  stage  in  place  of 
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TABLE  II 


Characteristic  length  divided  by  thickness 
divided  by  thickness  for  several  values  of 
divided  by  4ttM  . 

s 


vs.  strip  domain 
anisotropy  field 


V4™s  - 

00 

v» 

l/h 

0.50 

.0085 

0.60 

.0122 

0.70 

.0166 

0.80 

.0216 

0.90 

.0273 

1.00 

.0335 

1.50 

.0706 

2.00 

.1126 

3.00 

.1945 

4.00 

.2660 

5.00 

.3268 

6.00 

.3791 

7.00 

.4245 

8.00 

.4647 

9.00 

.5005 

10.00 

.5329 

15.00 

.6591 

20.00 

.7496 

4  2 


l/h 

l/h 

.0080 

.0076 

.0115 

.0110 

.0157 

.0149 

.0205 

.0195 

.0258 

.0246 

.0318 

.0304 

.0683 

.0660 

.1109 

.1088 

.1961 

.1967 

.2708 

.2746 

.3342 

.3406 

.3882 

.3965 

.4350 

.4446 

.4760 

.4866 

.5126 

.5239 

,5454 

.5572 

,6729 

.6861 

,7640 

.7776 
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the  sample.  The  resolution  of  the  microscope  (+0.5  to  +  0.2pm 
as  discussed  above)  is  generally  the  limitation  on  this 
calibration  so  that  1%  accuracy  ir.  a  measurement  of  100pm  is 
quite  possible. 

Experimentally,  the  domain  period,  P  ,  at  zero  field  can 

o 

be  measured  with  much  greater  accuracy  than  can  an  individual 
domain  diameter.  This  is  the  result  of  three  factors:  1)  The 
strip  domain  period  is  greater  than  the  bubble  diameters  by 
approximately  2  for  the  operating  bubble  size  and  4  for  the 
collapse  size.  2)  It  is  usually  possible  to  arrange  the  zero 
field  strips  into  a  parallel  array,  either  by  oscillating  them 
with  an  AC  bias  field  until  a  useful  pattern  develops  (this 
procedure  is  aided  by  the  application  of  an  in-plane  field 
component) ,  or,  as  suggested  by  Fowlis  and  Copeland,  by 
1  combing"  them  with  a  magnetic  probe.  Once  such  an  array  has 
been  formed,  the  distance  for  twenty  or  more  periods  can  be 
measured  and  an  average  value  computed  very  conveniently. 

3)  Effects  of  finite  coercivity  are  more  readily  minimized 
in  the  strip  pattern  than  for  bubble  domains  since  a  finite 
AC  bias  field  can  be  applied  and  gradually  brought  to  zero 
without  fear  of  collapsing  the  domains  of  interest.  By  the 
judicious  use  of  these  techniques  an  overall  accuracy  of  4 % 
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in  i  is  readily  achieved.  At  present  we  recommend  the  strip 
domain  measurement  at  H=0  as  the  most  convenient  and  accurate 
approach  to  the  determination  of  £. 
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9.  SATURATION  MAGNETIZATION,  4ttM  ,  AND  MAGNETIC  FIELDS 
_ s _ 

The  study  of  magnetization  probably  dates  back  to  the 
ancient  Greeks  or  Chinese  (Ref.  69) .  Since  that  time  a 
great  many  techniques  have  been  developed  for  its  measure¬ 
ment  and  many  of  these  have  been  applied  to  the  investigation 
of  orthoferrites  and  garnets.  As  the  development  of 
magnetic  bubble  domain  devices  has  progressed  toward  smaller 
bubbles  and  thinner  samples  some  of  these  techniques  have 
lost  their  usefulness,  particularly  for  thin  films  on  more 
massive  paramagnetic  substrates.  Others,  more  compatible 
with  the  magneto-optical  approach  used  for  other  parameters , 
have  been  developed.  In  the  present  report  we  will  touch 
on  these  earlier  methods  only  briefly,  saving  most  of 
our  attention  for  those  techniques  of  particular  use  in 
the  magnetic  bubble  area.  Furthermore,  excellent  reviews 
of  these  earlier  methods  already  exist  in  the  literature 
(Ref.  70,  71,  and  72). 

9 . 1  Traditional  Methods 

The  force  on  a  magnetized  sample  placed  in  a  non- 
uniform  magnetic  field  forms  the  basis  of  one  class  of 
magnetization  measurements  (Ref. 73  and  74).  Bozorth,  et  al 
(Ref.  75)  adapted  this  technique  to  the  low  temperature 
study  of  a  number  of  rare-earth  orthoferrites.  This  system 
involved  the  force  on  the  sample  in  the  field  of  an  electro¬ 
magnet  with  specially  formed  pole  pieces.  These  provided 
a  field  gradient  normal  to  the  net  field  direction  and 
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caused  a  force  on  the  sample  toward  the  higher  field  region. 
This  force  was  balanced  by  passing  a  current  through  a  small 
coil  riding  with  the  sample  and  just  balancing  the  sample 
magnetic  moment.  The  system  vras  calibrated  by  replacing  the 
sample  with  one  of  known  magnetization,  in  this  case  iron. 
Such  a  means  of  calibration  has  been  frequently  used. 

Force  systems  of  similar  types  have  been  employed  in  the 
study  of  garnets,  as  well  as  the  orthoferrites  (Ref.  76, 

7  7 ,  and  78). 

The  induction  of  current  in  a  coil  when  a  magnetic 
sample  is  moved  in  its  vicinity  dates  back  to  the  develop¬ 
ment  of  Faraday's  law  in  1834.  The  effect  has  been  used 
numerous  times  in  the  study  of  materials  presently  of 
interest  for  magnetic  bubble  applications  (Ref.  79  ,  76  , 

77,  and  80).  The  availability  of  lock-in  amplifiers  and 
other  adaptable  AC  electronics  has  caused  recent  interest 
toward  continuous  generation  of  an  alternating  induction 
signal.  In  these  either  the  sample  (Ref. 81  and  82)  or  the 
pick  up  coil  (Ref. 83  and  84)  are  vibrated  and  the  signal 
detected  by  means  of  a  lock-in  amplifier.  The  vibrating 
sample  system  has  been  used  by  Heinz,  et  al  (Ref.  85)  for 
garnet  films  on  paramagnetic  garnet  substrates  but, 
because  of  the  low  saturation  magnetization  of  practical 
magnetic  bubble  films,  accuracy  in  these  cases  is  limited. 

A  magnetized  sample,  when  placed  in  a  magnetic  field 
which  is  not  parallel  to  the  easy  magnetization  direction 
of  the  sample,  will  experience  a  torque.  In  applying  this 
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principle  to  the  measurement  of  magnetization  and  anisotropy, 

a  counter- torque  is  generally  applied  by  means  of  current 

carrying  coils  to  the  same  shaft  on  which  the  sample  is 

mounted.  Such  null  deflection  systems  with  automatic 

recording  features  have  been  described  in  detail  by 

Croft,  et  al  (Ref.  86),  by  Penoyer  (Ref.  87),  and  by  Byrnes 

(Ref.  88) .  An  instructive  schematic  diagram,  presented 

by  Chikazumi  (Ref.  89)  is  reproduced  as  Figure  14.  The 

paired  photocells,  B,  generate  an  unbalance  signal  when 

a  torque  is  applied  which  is  amplified  and  applied  to  the 

coil,  C,  in  such  a  way  as  to  restore  the  position  of  the 

entire  specimen,  coil,  mirror  assembly  to  very  nearly  its 

original  position.  The  current  required  is  a  measure  of 

the  torque  applied  and  is  plotted  vs .  the  direction  of  the 

applied  field  by  the  X-Y  recorder.  For  small  applied 

fields  (smaller  than  that  required  to  rotate  the  natural 

M  within  the  sample  significantly)  this  torque  is  a 
s 

measure  of  M  .  At  fields  high  enough  to  force  M  into 
s  s 

parallelism  with  H  the  torque  becomes  a  measure  of  the 
anisotropy  in  the  magnetic  energy  (see  Section  13) . 

The  torque  magnetometer  approach  has  been  used 
frequently  both  in  the  study  of  orthoferrites  (Ref.  90, 

91,  92,  and  93)  and  the  garnets  (Ref.  94,  95,  96,  97, 

and  98) .  Because  the  torque  method  is  applicable  to  the 
determination  of  magnetic  anisotropy  as  well  as  magnetization, 
several  of  these  references  as  well  as  others  cited  in 
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Section  17  deal  primarily  with  anisotropy.  Lacey  has 
recently  inaugurated  a  torque  system  in  which  the  effects 
due  to  paramagnetic  substrate  and  magnetic  film  can  be 
separately  determined  for  epitaxial  samples  (Ref. 99  and  100). 
There  are  sensitivity  problems,  however,  for  the  low 
anisotropy  materials  most  desirable  as  magnetic  bubble 
films  and,  in  any  case,  the  separation  of  film  and  substrate 
effects  requires  considerable  analysis  of  the  torque 
versus  applied  field  data. 

9.2  Magneto-optical  Methods 

Magneto-optical  approaches  to  the  determination  of 
magnetization  are  particularly  appropriate  for  thin, 
optically  transparent  bubble  domain  samples.  This  is 
especially  true  since  the  entire  fledgling  industry  depends 
upon  microscopic  examination  of  the  behavior  of  magnetic 
domains  via  the  Faraday  effect  for  the  evaluation  of 
both  materials  and  circuits  for  use  with  these  materials. 
Techniques  for  the  determination  of  magnetization  from 
optically  observed  behavior  in  various  magnetic  fields 
derive  from  the  same  theoretical  approaches  which  were 
central  to  the  determination  of  characteristic  length 
(Section  8) .  Either  the  behavior  of  magnetic  bubble 
domains  or  strip  domains  may  be  employed.  We  consider 
first  the  use  of  bubble  domains. 

The  theory  of  Thiele  (Ref. 29  and  30)  yields  the 
equilibrium  equation  for  bubble  domains  quoted  in  Section  8 
(Equation  8-1) •  Once  the  characteristic  length  and  thickness 
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have  been  determined  by  the  methods  discussed  previously  each 

pair  of  values  for  applied  field  and  bubble  diameter  serve 

to  determine  4nM  through  this  equation.  Hagedorn,  et  al 

(Ref.  58)  have  carried  out  a  least  squares  fit  to  this 

equation  of  careful  measurements  of  domain  diameter  and 

applied  field  (10  each)  and  conclude  that  relative 

determinations  of  47rMg  to  within  +  3%  can  be  made  using 

this  approach.  Note  that,  because  of  the  nature  of  this 

equation,  the  accuracy  of  4ttM  is  considerably  better 

s 

than  that  for  l  determined  in  the  same  way  (4  7%) .  It 

still  does  not  meet  the  requirements  laid  down  by  Bobeck, 

et  al  (Ref.  32)  of  1%  reproducibility  in  4ttM  for  devices. 

s 

It  also  has  all  of  the  difficulties  associated  with  it 
that  were  mentioned  in  Section  8  for  the  determination  of 
Si  from  bubble  dimensions. 

Fowlis  and  Copeland,  in  their  search  for  rapid  means 

of  evaluating  bubble  material  parameters,  have  suggested 

that  only  the  bubble  collapse  field,  Hq,  and  the  zero 

field  domain  period,  Pq,  be  measured  for  the  determination 

of  4ttMs  (Ref.  59).  They  quote  a  calculation  of  Goldstein, 

based  on  the  theories  of  Malek  and  Kambersky  (Ref.  61)  and 

Thiele,  which  relates  these  quantities  to  4ttM  .  The  results 

s 

are  shown  in  the  curve  of  Figure  15.  In  contrast  to  the 
technique  which  uses  only  bubble  measurements,  the  relative 
accuracy  of  the  4ttM  value  is  lover  than  that  for  J l  here, 
being  estimated  by  Fowlis  and  Copeland  as  7%  for  4rrM  and 


5%  for  l . 
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Fiqure  15.  Bubble  collapse  field,  H0,  divided  by  4ttMs  vs. 

strip  domain  period,  PQ,  at  H  *  0  divided  by 
thickness,  h.  After  Fowlis  and  Copeland,  ref.  59. 
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One  caution  concerning  the  use  of  the  bubble  collapse 
field  (and,  indeed,  all  bubble  measurements)  should  be  made 
here.  As  discussed  by  Bobeck  recently  (Ref. 101)  and 
observed  at  Monsanto  and  at  Hewlett  Packard  Companies 
(Ref. 102)  there  appears  to  be  more  than  one  type  of  bubble 
capable  of  at  least  metastable  existence  within  a  given 
bubble  film.  Bobeck  has  termed  all  of  those  which  are  not 
in  harmony  with  the  Thiele  theory  "hard"  bubbles,  since 
a  larger  magnetic  field  is  required  to  collapse  them. 

They  apparently  differ  from  "ordinary"  bubbles  primarily 
in  their  wall  energy  density.  Their  formation  seems  to  be 
favored  when  bubbles  are  formed  using  fast  rise-time  field 
pulses  such  as  are  frequently  employed  in  the  measurement 
of  mobility  (see  Section  11) .  Thus  one  can  be  misled  in 
the  measurement  of  collapse  field  appropriate  to  the 
Thiele  theory.  Generation  of  bubbles  by  methods  less 
energetic  than  fast,  high  field  pulses  are  to  be  preferred. 
Such  a  recommended  technique  would  be  the  cutting  of 
strip  domains  with  a  fine  magnetized  wire. 

One  further  magneto-optical  method  remains  to  be 
discussed,  namely,  the  use  of  the  full  hysteresis  curve 
(M/Mg  vs.  H) .  This  curve  is  readily  obtained  for  small 
regions  of  the  samples  by  mounting  a  photomultiplier  or 
other  optical  detector  at  the  eyepiece  of  a  polarizing 
microscope  equipped  with  bias  field  coils.  Cycling  the 
bias  field  to  saturation  while  recording  the  bias  coil 
current  versus  photomultiplier  output  using  an  X-Y 
recorder  yields  a  curve  like  Figure  16 .  A  quantitative 
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Hysteresis  curve  (M/Ms  vs.  H)  for  film  of  Gd 


M/M  vs.  H  curve  then  requires  only  a  knowledge  of  the 
bias  coil  calibration  and  reasonable  assurance  that  the 
M/M  curve  (for  which  M/M  =  1  is  usually  apparent)  is 
linear.  This  is  a  good  approximation  as  long  as  the 
polarizer-analyzer  pair  is  not  too  close  to  extinction. 

A  somewhat  more  elaborate  temperature  controlled  Faraday 
hysteresigraph  has  been  described  by  Beck  (Ref. 103)  and 
applications  of  a  similar  system  to  garnets  described  by 
Mee  (Ref . 104andl05) .  Both  of  these  systems  use  He-Ne  laser 
sources  while  we  have  found  that,  if  a  properly  equipped 
microscope  is  available  so  that  the  highly  collimated 
beam  is  not  required,  better  signal  to  noise  can  be  obtained 
using  an  incandescent  (quartz-iodine)  source  with  a  photo¬ 
multiplier  with  cutoff  in  the  red  (such  as  the  R.C.A.  6199 
with  S-ll  surface) . 

The  interpretation  of  the  hysteresis  curve  in  terms 

of  4ttMs  involves  recourse  to  theory,  such  as  that  of 

Kooy  and  Enz  (Ref.  62).  Craik  (Ref. 106)  has  outlined  a 

method  based  on  the  initial  slope  of  the  curve  and  unpublished 

curves.  Shaw,  et  al  (Ref.  60)  have  solved  Equations  9  and  10 

of  Kooy  and  Enz  for  the  values  of  H/4ttM  corresponding  to 

s 

0.1  intervals  in  M/M  .  The  results  for  high  anisotropy 
materials  (H  >>  4ttm  )  are  shown  in  Figure  17.  The 

3  S 

procedure  for  determining  M  then  involves  reading  values 

b 

of  M/M  and  H  from  the  hysteresigraph  for  a  sample  of  known 
s 

PQ  and  h.  The  origin  for  H  is  moved  to  the  right  on  the 
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Figure  17 


The  applied  field  H  necessary  to  achieve  various 

ze?oef^l^^rrati°S  “/Ms  f°r  samples  exhibiting 
field  strip  domain  repeat  distance  Pn. 

the  Les?Y  a"  EnZ'  N°te  the  normalizetions  of 
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hysteresis  loop  by  one  half  the  loop  width  at  H  =  0  and 

the  M/Mg  read  from  the  rising  portion  of  the  curve.  The 

theoretical  curves  then  permit  determination  of  H/4ttM  , 

s 

and  from  that  4irM  ,  for  a  particular  sample  at  a  particular 

s 

M/M  value.  Figure  16  illustrates  the  agreement  between 
s 

the  theory  (points)  and  an  experimental  curve  (continuous 
line) .  The  use  of  all  18  M/M  and  H  values  results  in  a 

o 

reasonably  accurate  value  of  4ttM  .  We  estimate  that 

relative  measurements  of  4tMs  can  be  made  from  sample  to 

sample  of  the  same  composition  with  this  technique  to 

within  +  3%.  The  term  "relative"  here  means  we  ignore 

uncertainties  in  index  of  refraction  (which  enters  the 

determination  of  h) ,  in  bias  coil  calibration,  in  graduated 

eyepiece  calibration,  and  in  approximations  in  the  theory 

since  these  errors,  if  present,  will  enter  each  determination 

of  4ttM  in  the  same  way. 
s 

Calculations  are  presently  in  progress  at  Monsanto 

(Ref.  68)  which  begin  from  the  more  accurate  Equations  1-3 

of  Kooy  and  Enz  and  take  into  account  finite  anisotropy 

down  to  Hn  =  2(4ttM  ).  Preliminary  results  indicate  that 

the  largest  changes  in  Figure  17  will  occur  for  large  M/M 

and  small  PQ/h  where  they  will  amount  to  an  increase  in 

H/4ttM  of  approximately  4%.  In  the  other  extreme  of  this 
s 

figure  the  changes  are  less  than  1%.  We  feel  that  the 
hysteresis  curve  approach,  utilizing  either  the  procedure 
just  outlined  or  a  direct,  least-squares  fit  to  the 
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means  for 


theoretical  equations,  is  the  most  accurate 
determining  the  magnetization  of  the  thin  epitaxial  films 
needed  for  magnetic  bubble  technology.  The  greater  convenience 
of  the  method  of  Fowlis  and  Copeland  will,  however,  cause 
It  to  be  used  for  more  routine  measurements. 
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10  DOMAIN  WALL  ENERGY 

The  two  previous  sections  outline  methods  for  the 

2 

determination  of  £  i  av/4lTMs  and  4itMs  •  Together,  then,  they 
provide  the  data  needed  to  calculate  the  domain  wall  energy 
per  unit  area,  o  .  Our  primary  concern  in  this  section  will 
be  the  discussion  of  other  methods  which  have  been  devised 
to  more  or  less  directly  determine  o  .  First,  however,  we 
briefly  mention  its  relationship  to  other  magnetic  parameters. 

The  form  of  domain  wall  which  normally  occurs  in  the 
garnets  (and  for  some  orientations  in  the  orthoferrites)  was 
first  investigated  by  Bloch  (Ref . 107andl08) .  In  this  Bloch 
wall  the  individual  spins  are  progressively  rotated  as  one 
traverses  from  one  domain  to  the  next,  the  spin  orientation 
always  lying  in  the  pl*"e  of  the  wall.  In  a  crystalline 
film  of  uniaxial  anisotropy  with  easy  axis  normal  to  the  film 
(a  first  approximation  to  the  situation  in  a  magnetic  bubble 
film)  the  net  magnetization  then  undergoes  a  180°  rotation 
between  adjacent  domains.  Two  terms  contribute  to  the  increase 
in  energy  which  this  configuration  causes,  being  offset  by 
the  decrease  in  magneto-static  energy  resulting  from  the 
breakup  into  domains.  These  two  terms  arise  from  isotropic 
exchange  and  uniaxial  anisotropy  characterized  by  the  constants 
A  and  Ku,  respectively.  Upon  minimizing  the  total  energy 
with  respect  to  the  rate  of  rotation  of  the  local  moment,  one 
finds  a  total  wall  energy  per  unit  surface  area  given  by  a 
combination  of  isotropic  exchange  and  anisotropy  as 
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(10-1) 


°w  =  4 (AKu» 1/2 

The  width  of  the  wall  (defined  as  that  width  required  to 
encompass  the  total  180°  rotation  if  the  maximum  rate  of 
rotation  at  the  wall  center  were  maintained  throughout)  is 
given  in  this  same  context  as 

#.  =  tt  ( A/K  ) 1/2  (10-2) 

w  u 

As  already  mentioned  the  anisotropy  constant  is  related  to 
the  field  required  to  force  the  magnetization  into  the  hard 
direction,  H,,  as  K  =  H,M  /2.  Thus  the  present  section, 
together  with  those  dealing  with  M  and  H  (Sections  9  and 

S  A 

13)  permit  the  complete  characterization  of  magnetic 
bubble  films  in  terms  of  these  basic  parameters. 

Kurtzig  and  Shockley  (Ref .109)  developed  a  method  for 
the  determination  of  domain  wall  energy  in  the  orthoferrites 
which  involves  the  increase  in  wall  length  brought  on  by  the 
application  of  a  particular  type  of  inhomogeneous  field. 

They  considered  straight  wall  strip  domains  with  an  array  of 
parallel  wires  laid  across  the  domains  at  right  angles. 
Application  of  equal  but  opposite  current  to  adiacent  wires 
of  this  set  caused  a  spatially  periodic  distortion  of  these 
straight  walls.  Although  the  calculations  and  experiments 
were  carried  out  under  conditions  permitting  several 
simplifying  assumptions  (the  conductors  were  far  enough 
from  the  sample  to  yield  only  spatially  sinusoidal  field 
variations,  the  amplitude  of  wall  motion  was  small  compared 
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to  width,  the  domain  walls  always  contained  the  easy  axis  of 

magnetization)  the  expressions  remained  formidable  and  will 

not  be  reproduced  here.  The  method,  averaged  over  ten  deter- 

2 

minations  on  five  samples  of  ErFeO^,  yielded  1.95  ergs/cm  + 
5%.  The  tota]  variation  among  the  ten  cases  cited  was  14%. 
For  the  materials  with  smaller  domain  dimensions  presently 
of  interest  this  approach  would  be  considerably  less  accurate 
since  much  smaller  domain  motion  would  need  to  be  measured. 

Nemchik  and  Charap  (Ref . 110andlll)have  investigated  the 
behavior  of  an  array  of  bubble  domains  in  GdIG  and  developed 
two  additional  means  of  obtaining  o^.  In  this  case  a 
sufficient  density  of  bubbles  was  formed  by  means  of  a  pulse 
coil  to  inhibit  the  tendency  to  run  out  into  strip  domains 
at  low  fields.  Thus,  even  at  zero  bias  field,  the  domain 
array  was  that  of  nearly  circular  bubbles.*  Charap  and 
Nemchik  have  also  made  use  of  the  fact  that  GdIG  has  a 
compensation  temperature  near  room  temperature  (Tc  ■=  9.6°C) 
(Ref. 105).  This  leads  to  a  large  variation  with  temperature 
of  nearest  neighbor  bubble  spacing  and  bubble  size  at  zero 
applied  bias  field.  By  fitting  variable  temperature  data  to 
theoretical  curves  and  assuming  that  the  domain  wall  energy 


*  Cape  and  Lehman  (Ref  H2andll3)  and  Kaczer  (Ref .  3 14)  have 
considered  the  energetics  of  bubble  and  strip  domain 
patterns.  They  showed  that  the  energy  difference  between 
such  a  bubble  array  and  the  strip  pattern  is  not  large 
and  discussed  the  barriers  separating  the  two  patterns. 
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itself  is  not  a  function  of  temperature  near  T  ,  they  have 

c 

derived  values  of  from  each  of  these  temperature 
variations.  Cape  and  Lehman  point  out  an  error  in  their 
computation,  however,  so  that  at  the  present  time  we  can 
only  say  that,  done  properly,  the  method  appears  capable 
of  ^  20%  accuracy  in  materials  having  significant  tempera¬ 
ture  coefficients  in  the  region  of  interest.  Such 
variability  with  temperature  is  not  desirable  in  a  device 
film,  of  course. 

In  considering  the  alternative  approaches  to  the 

measurement  of  domain  wall  energy  per  unit  area,  the  best 

balance  of  speed,  convenience  and  accuracy  is  found  in  the 

determination  of  £  and  4irMs  by  the  methods  already  discussed 

2 

and  the  calculation  of  a  -  £4ttM  . 

w  s 


11.  DOMAIN  WALL  DYNAMIC  PROPERTIES 


The  parameter  most  commonly  used  today  for  the  specifi¬ 
cation  of  dynamic  domain  wall  response  is  the  domain  wall 
mobility,  yw«  This  is  defined  in  terms  of  the  wall  velocity, 
v,  the  field  required  to  first  cause  wall  motion,  Ho,  and 
the  instantaneous  applied  field,  H,  as 

v 


yw  H-H 


(11-1) 


o 

This  parameter  has  direct  bearing  on  the  speed  with  which 
bubbles  can  be  moved  in  a  device  and  therefore  on  the  maximum 
data  rate  possible.  As  work  has  progressed,  however,  it 
has  become  clear  that,  in  high  mobility  materials,  a  strict 
proportionality  between  domain  wall  velocity  and  net  drive  field 

holds  only  for  small  H-H  (Ref .115andll6 ) .  For  this  reason 

o 

as  well  as  others  to  be  mentioned  later,  the  present  section 
will  extend  beyond  the  measurement  of  mobility  only. 

Attempts  to  relate  mobility  to  other  parameters 
characteristic  of  the  materials  under  study  have  met  with 
limited  success.  An  excellent  review  of  this  situation  has 
recently  been  presented  by  Hagedorn  (Ref.  63)  and  it  is  not 
the  function  of  the  present  publication  to  repeat  these 
details.  The  reader  is  also  referred  to  the  earlier  reviews 
of  Dillon  (Ref. 117)  and  Kittel  and  Galt  (Ref .1184*  The  one 
relationship  which  appears  to  be  well  established  is  a  pro¬ 
portionality  of  mobility  and  domain  wall  width  as  given  by 
equation  (10-2) . 


1J--1  Methods  for  Bulk  Samples 


The  classic  method  (since  about  1950)  for  the  measurement 
of  domain  wall  motion  has  been  discussed  in  considerable 
detail  by  Galt  (Ref. 119).  It  involves  cutting  a  "picture- 
frame"  ring  with  each  leg  of  the  sample  oriented  along  a 
direction  of  easy  magnetization.  Such  a  structure  is  illus¬ 
trated  in  Fig.  18.  in  this  case  the  domain  wall  can  be 
propagated  through  the  sample  in  a  direction  normal  to  the 
plane  of  the  "picture-frame"  by  means  of  a  coil  wound  on  one 
leg  of  the  frame.  The  motion  of  the  domain  wall  can  be  detec¬ 
ted  by  the  emf  generated  in  another  similar  coil.  Considerable 
care  must  be  exercised  to  insure  that  only  one  wall  of  the 
moving  type  is  present.  This  sometimes  requires  careful 
proportioning  of  the  sample  in  order  to  favor  a  particular 
wall  (Ref. 120).  The  wall  velocity  is  then  determined  by 
applying  a  pulse  of  current  to  the  drive  coil  and  measuring 
the  magnitude,  V,  of  the  voltage  pulse  generated  in  the 
pick-up  coil.  Ideally  this  is  a  flat,  square  pulse  from 
which  the  wall  velocity  follows  using  Faraday's  law  and  the 
system  geometry: 


v 


_ V 

N8ttM  W  x  10 
s 


(11-2) 


Here  N  is  the  number  of  turns  in  the  pick  up  coil,  W  is  the 
dimension  of  the  sample  leg  perpendicular  both  to  its  axis 
and  the  direction  of  wall  motion,  and  the  units  are  mixed: 
v-cm/aec;  V-volts,  Mg-Gauss,  and  W- cm.  The  mobility  then 
follows  from  equation  (11-1) . 
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Tms  technique,  without  the  use  of  a  full  "picture- 
frame"  sample,  has  been  applied  to  yttrium  orthoferrite  by 
Umebayashi  and  Ishikawa  (Ref. 121)  for  temperatures  down  to 
77  K.  Yttrium  iron  garnet  was  measured  using  the  full 
frame  configuration  by  Hagedorn  and  Gyorgy  (Ref. 122,  and  by 
Wanes  (Ref. 123)  down  to  4.2»K.  A  difference  of  a  factor  of 
4-7  between  these  results  probably  resulted  from  the 


presence  of  multiple  domain  walls  in  the  earlier  measurements, 
wanes  also  observed  the  decrease  in  mobility  brought  about 
by  the  addition  of  small  amounts  of  Si,  a  clear  indication 
to  today's  crystal  grower  of  the  need  for  purity. 


+  ReSponse  of  sin^  crystal  toroids  of  YIG  doped  with 
Yb  and  Er3+  ions  has  been  investigated  by  Harper  and 

Teale  (Ref. 124  andl25, .  Because  of  the  sample  geometry  used 
they  suspect  the  presence  of  multiple  domain  walls  and  there¬ 
fore  work  with  switching  times  rather  than  mobilities.  The 
temperature  dependences  observed  (large  peaks  in  switching 
time  below  approximately  100-K)  are  certainly  reflections  of 

the  damping  introduced  by  the  Yb  and  Er  rare  earth  ions, 
however. 


While  this  "picture- frame”  technique  is  an  excellent 
one  for  the  study  of  domain  wall  dynamics  it-,  chief  drawback, 
sample  preparation,  should  be  obvious.  Such  complex  samples 
cannot  be  formed  from  the  thin  films  needed  for  devices  nor, 
in  fact,  could  they  be  formed  at  all  from  uniaxial  materia/ 
and  still  maintain  the  easy  axis  orientation  in  all  legs. 


Thus  the  approach  may  remain  quite  useful  for  fundamental 
studies  but  has  no  place  in  the  routine  characterization  of 
films  for  bubble  devices. 

It  is  not  necessary  to  have  a  closed  magnetic  loop  such 
as  the  "picture-frame"  geometry,  of  course,  in  order  to  induce 
an  electrical  signal  as  a  result  of  domain  wall  motion.  Asti, 
et  al  (Ref. 126, 127  andl28)  have  made  use  of  ellipsoidal 
samples  placed  within  a  coil  which  is  in  turn  within  a  larger 
solenoid.  A  pulsed  field  is  generated  by  means  of  a  capaci¬ 
tor  discharge  through  the  solenoid  and  a  signal  proportional  to 
the  rate  of  change  of  sample  magnetization  is  developed  in 
the  pick  up  coil.  Effects  not  associated  with  the  sample 
are  bucked  out  by  a  matched  compensating  coil  also  within 
the  solenoid.  A  number  of  ferrites  have  been  studied  in 
this  way  but,  once  again,  application  to  thin  films  is  not 
convenient. 

Vella-Coliero,  et  al  (Ref. 129  and!30)  have  measured  the 
real  and  imaginary  parts  of  the  susceptiblity  of  small,  bulk 
crystals  of  garnets  using  a  modification  of  a  nuclear  magne¬ 
tic  resonance  system  (Ref .131).  They  determine— a  relaxation 
frequence,  ooc,  at  which  the  imaginary  term  in  the  susceptibi¬ 
lity  has  a  maximum  while  the  real  part  has  fallen  to  1/2  of 
the  low  frequency  limit,  xo*  The  domain  wall  mobility  is 
then  calculated  as 
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where  d  is  the  domain  width  in  the  sample.  For  two  series 
of  Tb  and  Dy  additions  to  Gd3Fe5C>12  agreement  to  within  10% 
with  the  method  of  Bobeck  (see  below)  is  indicated  but  no 
indication  of  absolute  accuracy  is  given. 

In  the  second  of  these  papers  (Ref. 130)  the  authors 
have  investigated  a  series  of  high  mobility  Y,  Gd  Fe^O. 0 
samples  and  have  observed  a  resonant  behavior  near  100  MHz. 
Such  a  resonance  implies  a  domain  wall  inertia  and  thus  an 
effective  mass  for  the  wall.  The  mass  per  unit  area  was 
calculated  from  the  resonant  frequency  w  as 


4M 
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(11-4) 


The  results  indicate  that  m  rises  from  ^  2  x  10*11  to 
-11  2 

12  x  10  gms/cm  as  the  yttrium  content  increases  and 
clearly  demonstrates  a  means  of  evaluating  this  wall  mass  in 
suitable  materials.  In  addition,  mobility  values  calculated 
from  equation  (11-3)  increase  from  'v  5000  to  12000  cm/sec  Oe 
for  the  same  composition  change.  These  are  the  highest  room 
temperature  mobilities  reported  to  date.  Henry  (Ref. 132) 
has  pointed  out  that  wall  mass  effects  may  become  significant 
in  device  applications  of  such  materials. 


11-2  Pulsed  Field  Methods  for  Thin  Film  or  Platelet  Samples 
The  application  of  a  pulsed  field  to  the  thin  platelet 
or  film  samples  of  interest  for  bubble  devices  does  lead  to 
measurable  effects.  Middelhoek  (Ref. 133)  has  described  a 
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system  appropriate  for  Ni-Fe  films  in  which  wall  position 
was  detected  via  the  Kerr  effect.  The  change  in  position 
of  a  single  wall  in  an  initially  saturated  film  following 
a  series  of  field  pulses  was  interpreted  in  terms  of  a 
velocity.  The  mobility  followed  as  the  slope  of  the 
velocity  ve.  pulsed  field  amplitude  curve.  Other  results  in 
Ni-Fe  films  have  recently  been  reported  by  Konishi  (Ref.134). 
In  magnetic  bubble  materials,  domain  walls  are  usually 
detected  via  the  Faraday  effect.  Seitchik,  et  al  (Ref. 135) 
have  discussed  a  technique  used  in  the  orthoferrites  in 
which  the  variation  of  light  passing  through  the  sample 
(and  passing  through  the  usual  polarizer-analyzer  pair)  is 
measured  as  a  short  risetime  field  pulse  is  applied.  In 
order  to  minimize  inductive  effects  a  small  (20  turns,  1.8mm 
diameter)  coil  was  used  and  sampling  and  time  averaging 
techniques  applied  to  the  resulting  signal  to  improve  the 
signal  to  noise  ratio.  The  relaxation  time,  t,  of  the 
domain  walls  to  their  new  configuration  in  the  pulsed  field 
must  be  considerably  longer  than  the  rise  time  of  the  applied 
pulse  and  the  field  pulse  amplitude  must  be  less  than  1/2  of 
the  saturation  field  to  maintain  accuracy.  In  this  case 
starting  with  strip  domains 

yw  =  <VHo>  t_1  dl-5) 

where  XQ  is  the  equilibrium  displacement  of  a  wall  in  an 
applied  field  Hq.  Shaw,  et  al  (Ref.  39)  have  discussed  a 
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means  of  calculating  Xq/Hq  from  a  knowledge  of  the  h,  p  , 
and  4trMs  using  Fig.  17.  At  Monsanto,  the  Seitchik  method 
has  been  used  extensively  in  mobility  measurements  in  thin 
film  magnetic  garnets.  Efforts  to  extend  this 
method  to  high  mobility  materials  have  yielded  values  to 
approximately  l,000cm/sec  Oe  in  garnets  from  2  turn  or  5 
turn  coils  which  allow  pulse  risetimes  of  1  and  2  nsec, 
respectively.  A  mercury-wetted  reed  pulser  (Spencer  Kennedy 
Laboratories  Model  503A)  with  0.5  nsec  risetime  is  used  to 
drive  these  coils  in  a  strip-line  configuration.  Accuracy 
is  estimated  at  +  25%  in  the  high  mobility  range  while  the 
value  of  the  mobility  is  a  function  of  pulse  amplitude  as 
expected  from  the  results  of  Calhoun,  et  al  (Ref. 115)  and 
Bonner,  et  al  (Ref. 116).  Shumate  (Ref.  47  and  40)  and  Shaw,  et  al 
(Ref.  39)  have  discussed  the  use  of  this  approach  in  the 
garnets.  Nemchik  and  Charap  (Ref .136)  have  discussed  a 
similar  approach  which,  however,  starts  with  a  zero  field 
bubble  array  rather  than  strip  domains. 

Another  approach  to  the  mobility  which  also  employs  a 
pulsed  field  and  Faraday  rotation  is  the  bubble  collapse 
method  of  Bobeck  (Ref.137).  More  detailed  analysis  of  the 
method  (but  including  some  approximations)  has  been  presented 
by  Callen  and  Josephs  (Ref.138).  The  technique  involves  the 
application  first  of  a  DC  bias  field  which  stabilizes 
isolated  bubble  domains.  One  of  these  is  now  subjected  to 
a  further  bias  pulse  driving  it  toward  collapse.  Although 
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this  total  field  (Hb^as  +  Hpulse^  may  be  ^ar9e  enough  to 
cause  collapse  if  applied  as  a  constant  bias  field  for  a 
pulse  of  sufficiently  short  duration  the  bubble  will  still 
be  present  and  return  to  its  original  size  after  the  pulse 
terminates.  This  is  a  result  of  the  finite  speed  with 
which  the  domain  wall  can  move  toward  collapse  during  the 
pulse.  The  pulse  amplitude  and  duration  at  which  collapse 
just  takes  place  (i.e.  reducing  either  quantity  results  in 
the  continued  existence  of  the  bubble)  forms  the  basis  of 
the  measurement.  It  is  important  to  realize  that  the 
minimum  instantaneous  diameter  from  which  a  bubble  can 
recover  is  smaller  than  the  static  diameter  for  collapse, 
being  approximately  one  half  of  that  static  dimension  for 
a  bias  field  equal  to  the  strip  to  bubble  transition  field. 
The  original  Thiele  paper  (Ref.  29)  discusses  this  dynamic 
stability  limit  but  more  explicit  consideration  of  the 
limit  as  it  applies  to  the  present  method  has  been  given 
elsewhere  (Ref. 139  andl40) .  Shumate  (Ref.  47)  has  reported 
Rossol's  observations  of  the  dynamically  allowed 
range  of  diameters  using  stroboscopic  techniques.  If  the 
pulse  is  applied  to  a  bubble  initially  at  the  strip  to 
bubble  transition  the  total  distance  traveled  by  a  wall  to 
just  reach  collapse  is  Ar  =  rg_B  -  ro/2.  Bobeck  suggests 
plotting  the  inverse  of  the  pulse  duration,  T,  necessary  to 
bring  about  bubble  collapse  vs^.  the  field  pulse  amplitude 
Some  of  his  original  data  (Ref. 137)  is  reproduced  in  Fig.  19 
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Figure  19.  Bubble  collapse  data  for  a  uniaxial  garnet. 
After  Bobeck,  ref.  137. 


(11-6) 


The  mobility  is  then  calculated  as: 

_  (Ar)  A (1/T) 

“w  ”  iHp 

The  approach  of  Callen  and  Josephs  (Ref. 138)  is  couched 
in  terms  of  fields  and  pulse  lengths  only  (with  all  variables 
elegantly  normalized) .  For  pulse  fields  considerably  above 
those  required  to  cause  collapse  their  method  yields 

(V)'1  £  al  +  a2  Hp  (11-7) 

where  a^  and  a 2  are  functions  of  Hc  ,  4ttMs#  H  ,  and  h.  This 
can  be  cast  in  the  same  form  as  the  Bobeck  approach  by 
differentiation  in  which  case 


1_  A (1/T) 
a2  ~1h  * 


(11-8) 


In  the  present  notation 
-1 


2  (4ttMr)  h{[l-3/4  U/h)  +  (H1/4ttMs)] 
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which  replaces  Ar  in  the  Bobeck  expression.  Here  is  the 
DC  bias  field  applied. 

This  approach  has  been  widely  used  and  has  the  advantages 

of  yielding  results  on  a  wide  variety  of  materials  and 

requiring  a  minimum  of  additional  equipment.  Calhoun »  et  al 

(Ref. 115)  have  used  it  to  study  the  variation  of  velocity 

with  pulse  amplitude  in  high  mobility  Ga-YIG  where  they  find 

u  ^  fiT  .  Such  a  phenomenon  suggests  that  the  mobility 
v  p 

(usually  measured  at  small  pulse  amplitudes)  is  not  the  most 


93 


useful  parameter  for  predicting  device  performance.  Velocity 
measured  at  a  pulse  amplitude  characteristic  of  the  contem¬ 
plated  device  would  be  much  more  realistic  as  a  starting 
point  in  such  a  prediction.  Several  authors  (Ref.  141  &  142)  have 
recently  reported  results  on  various  Y-RE  garnets  indicating 
a  velocity  saturation  effect.  Both  the  bubble  collapse 
(Ref. 141)  and  step  function  response  (Ref. 142)  methods  were 
used.  Slonczewski  (Ref. 143)  has  attacked  this  problem  on 
a  theoretical  basis. 

There  are  a  number  of  areas  of  concern  in  the  application 
of  the  bubble  collapse  method.  First  (and  probably  least 
important)  the  dynamic  collapse  radius  mentioned  above  is  a 
function  of  the  DC  bias  field,  being  rc/2  only  at  the  strip 
to  bubble  transition  field  and  approaching  tq  as  the  DC  bias 
is  raised  toward  the  collapse  field.  Second,  Cape  (Ref.  144) 
has  pointed  out  that  the  Bobeck  expression  actually  involves  a 
linear  approximation  to  an  exponential  (the  exponential 
approach  of  the  bubble  radius  to  its  new  equilibrium)  and,  as 
such,  requires  pulse  fields  considerably  greater  than  the 
minimum  necessary  to  bring  about  collapse.  The  same  idea  is 
involved  in  producing  the  approximation,  equation  (11-8), 
from  the  theory  of  Callen  and  Josephs  above.  Cape  also  has 
some  more  subtle  restrictions  on  the  Bobeck  method  which 
would  seem  to  require  that  the  initial  radius  (and  therefore 
H^)  must  be  quite  close  to  the  collapse  value.  If  this  is 
true  the  accuracy  of  the  method  would  be  seriously  impaired. 


Third,  the  matter  of  different  kinds  of  bubbles  ("hard" 
bubbles)  raised  earlier  becomes  particularly  troublesome 
here.  The  bubble  under  study  must  be  annihilated  in  such  a 
measurement  and  it  is  very  tempting  to  reduce  the  bias  and 
pulse  the  coil  hard  in  order  to  generate  new  bubbles. 
Unfortunately  this  seems  to  be  the  very  best  way  to  generate 
these  anomalous  bubbles  which  do  not  behave  in  quantitatively 
the  same  way  that  "ordinary"  bubbles  do  and  lead  to  consider¬ 
able  scatter  in  the  data.  The  best  way  to  avoid  this  problem 
seems  to  be  each  time  to  cut  the  strips  with  a  fine  wire  as 
mentioned  earlier. 

11.3  Methods  Utilizing  the  Isolated,  Straight  Domain  Wall 

In  the  orthoferrites  an  isolated  wall  can  be  generated 
by  applying  a  sufficiently  strong  gradient.  Such  an 
arrangement  has  been  used  in  the  study  of  mobility  by  Rossol 
(Ref. 145  andl46)  and  Shumate  (Ref. 147).  Rossol's  approach 
has  been  to  apply  an  alternating  field  in  addition  to  the 
gradient  and  stroboscopically  determine  the  total  displace¬ 
ment,  X,  of  the  wall.  As  the  frequency  of  the  constant 
amplitude  alternating  field  is  raised  the  wall  velocity 
increases  until  it  becomes  limited  by  the  mobility  of  the 
material.  The  angular  frequency,  u)c,  at  which  the  wall 
displacement  has  fallen  to  1/  /2  of  its  low  frequency  value, 
Xo,  is  determined  from  a  frequency  vs.  displacement  plot. 

The  mobility  is  then  given  by 


Rossol  has  reported  mobility  for  a  number  of  the  rare-earth 
orthoferrites  including  values  as  high  as  b0,000cm/sec.  Oe 
in  carefully  prepared  YFeO^  at  77°K.* 

Shumate  (Ref. 147  andl49)  has  used  a  similar,  straight 
wall  technique  in  which  the  alternating  field  is  replaced 
by  a  square  pulsed  field  of  duration  such  that  the  domain 
wall  does  not  reach  its  new  equilibrium  before  the  pulse  is 
removed.  The  wall  approaches  its  new  position,  Xq,  during 
the  pulse  according  to  the  equation 

X(t)  =  XQ  [1  -  exp  (-yH  t/XQ)]  (11-11) 

Thus  measurements  of  the  maximum  displacement  for  various 
pulses  of  known  duration  and  amplitude  serve  to  determine 
u  .  Steady  rather  than  stroboscopic  illumination  has 

W 

generally  been  employed.  For  the  orthoferrites  this 
technique  appears  to  be  the  most  accurate  available,  the 
reported  reproducibility  on  a  given  sample  of  high  quality 
being  several  parts  per  thousand  and  the  absolute  accuracy, 

5%.  Shumate  investigated  the  anisotropy  of  mobility  - 
arising  from  the  presence  of  Bloch  and  Neel  walls  for  different 
wall  orientations  in  the  orthof errifes  (Ref. 150)  -  and  the 


*  Marsh,  et  al  (Ref. 148)  have  briefly  described  a  related 
technique  in  which  an  isolated  bubble  is  modulated  in 
size  by  an  AC  bias  component  and  the  field  amplitude 
required  for  constant  amplitude  in  radius  is  measured. 
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changes  of  mobility  brought  about  by  annealing  and  x-ray 
irradiation  (Ref. 149  andl51) .  (See  also  ref. 152) . 

The  difficulties  in  generating  the  field  gradient 
needed  to  yield  a  single  straight  wall  in  the  magnetic  bubble 
garnets  have  already  been  discussed  (see  Section  7) .  This 
limitation  severely  restricts  the  application  of  this 
technique  for  today's  magnetic  bubble  materials.  So  far  as 
the  authors  are  aware  no  analysis  of  the  corrugated  or  finger 
like  wall  which  occurs  in  the  garnets  at  reasonable  gradient 
fields  has  been  attempted  so  that  the  approach  cannot  be 
recommended  for  general  use. 

11.4  Bubble  Translation  Methods 

Several  schemes  have  been  proposed  for  the  study  of  the 
dynamic  translational  response  of  bubbles.  Such  a  measure¬ 
ment  is  to  be  desired  in  that  cylindrical  domain  walls  are, 
in  fact,  the  ultimate  moving  element  in  this  subject,  not 
individual  strip  domain  walls.  Heinz,  et  al  (Ref.  35)  (see 
also  ref. 153)  have  employed  conductor  loops  deposited 
directly  onto  the  bubble  film  surface.  These  were  pulsed 
so  that  a  net  field  difference  was  established  which  tended 
to  transfer  a  bubble  from  one  to  the  next.  Approximating 
this  as  a  constant  gradient  AH/Ax  between  the  loops  and 
knowing  the  minimum  pulse  duration,  tQ,  required  to  transfer 
a  bubble  of  diameter,  d,  between  loops  of  spacing  r  they 
were  able  to  calculate  the  approximate  domain  mobility 

=  _ r _ 

dtQ  (AH/Ax)  (11-12) 
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Here  „d  =  Pw/2.  Such  a  system  is  useful  in  the  garnets  but 
requires  the  deposition  of  the  conductor  pattern  (or  the 
close  approach  as  such  a  pattern  on  a  separate  microscope 
slide)  as  well  as  the  above  approximations  for  its  applica- 

tion. 

Copeland  and  Spiwak  (Ref. 154)  have  employed  a  single 
straight  conductor  to  generate  the  translating  gradient  field. 
They  present  theoretical  results  indicating  that  the  driving 
force  is  constant  to  within  20%  for  bubble  center  positions 
within  half  of  the  bubble  radius.  They  also  consider 
various  conductor  to  sample  spacings  and  sample  thicknesses. 
In  the  experiment  the  bubble  is  driven  using  a  bipolar  square 
wave  so  that  it  remains  within  the  region  mentioned  (two 
other  conductors  carrying  direct  current  are  also  present  to 
contain  the  bubble) .  In  arriving  at  mobility  values  these 
authors  also  identify  the  drive  field  HF  as  Hp  =  R  AH/AX 
for  a  bubble  of  radius  R.  Mobilities  obtained  by  this  method 
agree  to  within  10%  with  those  measured  by  Shumate  using  the 
straight  wall  method.  Unfortunately  the  method  appears  to  be 

limited  to  the  orthoferrites  (Ref. 155). 

Rossol  has  applied  the  stroboscopic  approach  to  the 

study  of  bubble  motion  under  a  permalloy  disk  (Ref. 156) 
and  in  a  simple  T-bar  circuit  (Ref. 157)  (again  on  ortho¬ 
ferrites).  In  the  former,  the  phase  lag  between  an  in-plane 
rotating  field  and  the  bubble  motion  was  measured  as  a  func¬ 
tion  of  the  rotation  frequency.  After  considerable  analysis 
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of  the  disk-bubble  interaction,  the  mobility  was  calculated 
and  agreed  well  with  that  found  by  the  straight  wall 
technique.  In  the  T-bar  circuit,  the  time  delay  between 
various  positions  for  a  circulating  bubble  in  the  device  was 
measured  and  a  velocity  obtained  for  each  region  of  the  device. 
In  this  case  the  field  gradients  and  corresponding  mobility 
again  required  detailed  analysis  which  was  not  presented.  The 
method  is,  in  fact,  of  more  usefulness  in  the  design  and 
analysis  of  magnetic  overlays  for  devices  than  for  material 
characterization. 

Vella-Coleiro  and  Tabor  (Ref. 158)  have  recently  proposed 
a  pulsed  gradient  system  of  considerable  accuracy  for  use  in 
the  garnets.  This  involves  the  used  of  two  parallel  conductors 
carrying  current  in  the  same  direction  which,  when  pulsed, 
yield  a  gradient  field  at  their  center  which  drives  a  bubble 
perpendicular  to  their  direction.  In  order  to  reduce  the 
effect  of  the  decreasing  field  seen  by  the  bubble  as  it  moves 
in  the  pulsed  field  they  superimpose  a  linearly  increasing 
field  component  by  means  of  an  additional  conductor  in  the 
form  of  a  rectangle.  This  ramp  contribution  to  the  field 
must  be  computed  and  properly  adjsuted  for  each  measurement 
of  bubble  motion  vs.  pulsed  field  gradient.  Further,  the 
total  travel  is  limited  to  ^  10%  of  the  distance  between 
the  conductors  in  order  to  remain  in  a  region  of  reasonably 
constant  gradient  3%).  For  the  circuit  described  this 
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total  travel  is  5um  so  that  its  measurement  will  pose  a 
significant  limitation  on  the  overall  accuracy  of  the  method. 
One  further  limitation  comes  from  the  decrease  in  gradient 
as  the  distance  from  the  plane  of  the  conductors  increases, 
being  approximately  13%  at  10y  and  limiting  the  technique  to 
samples  of  about  this  thickness.  Within  all  of  these 
limitations  the  method  yields  the  mobility  (and  the  dynamic 
coercivity,  Hc)  from  the  equation 

v  =  AX/At  =  1/2  yw  -  f  Hc^  (11-13) 

where  I  is  the  current  in  the  parallel  conductors  and  D  is 
their  spacing.  Several  measurements  under  various  pulse 
conditions  are  of  course  necessary  to  determine  both  nw  and 
Hc.  Cape  (Ref. 144)  has  proposed  a  similar  system  using  two 
wires  without  the  ramp  field  refinement.  He  has  also  dis¬ 
cussed  some  of  the  real  world  difficulties  which  occur  with 
such  a  system. 

11.5  Conclusions 

In  closing  this  section  on  the  many  approaches  which 
have  been  made  to  the  measurement  of  domain  wall  dynamics, 
several  points  should  be  made. 

1.  The  relationship  between  domain  wall  velocity  and 
drive  field  is  not  always  linear,  particularly  for 
materials  composed  of  ions  with  nearly  zero  orbital 
angular  momentum.  The  highest  mobilities  occur  for 
drive  fields  which  are  too  small  to  be  practical.  At 
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useful  drive  fields  the  domain  velocity  saturates  in 
many  cases  which  results  in  a  low  effective  mobility. 
Therefore  the  concept  of  a  mobility  loses  much  of  its 
usefulness.  If  mobility  values  are  quoted  they  should 
be  accompanied  by  a  drive  field  value  or  by  an 
assurance  that  the  value  is  independent  of  drive  field 
over  a  specified  range.  The  preferable  approach  would 
seem  to  be  to  quote  velocity  vs^.  drive  field  values  or 
curves . 

2.  It  would  appear  at  this  point  that  the  bubble 
translation  method  of  reference  158  is  the  best  approach 
to  the  accurate  measurement  of  velocity  vs.  drive 
field  while  the  step  function  response  of  linear  walls 
is  most  convenient  as  a  quality  control  test  in  a  pro¬ 
duction  situation.  In  spite  of  the  many  techniques  which 
have  been  proposed,  however,  it  seems  likely  that  the 
ideal  method  for  routine  evaluation  of  dynamic  properties 
of  bubble  materials  is  still  to  come. 

3.  Even  when  reliable  mobility  or  velocity  data  are  in 
hand,  the  prediction  of  device  operating  speeds  and 
margins  remains  a  formidable  and  imprecise  task.  Several 
useful  discussions  have  been  presented,  however,  and  are 
listed  here  as  ref erencesll>9 , 160,  andl61. 


101 


12.  FERROMAGNETIC  RESONANCE  (FMR) 


This  phenomenon  can  be  pictured  in  terms  of  the 
precession  of  the  total  magnetic  moment  of  the  specimen 
around  the  direction  of  an  applied  magnetic  field.  When 
a  transverse  field  is  applied  at  the  frequency  of  this 
precession,  energy  is  absorbed  as  the  system  is  excited 
(in  a  quantum  mechanical  picture)  to  higher  Zeeman  levels. 
YIG  and  the  related  rare-earth  garnets  present  a  nearly 
ideal  system  for  the  study  of  this  phenomenon  for  they 
are  good  insulators,  possess  strong  FMR  lines,  and  are 
theoretically  tractable.  Because  of  this  and  because 
commercial  microwave  applications  exist  for  YIG  the  litera¬ 
ture  on  this  subject  is  enormous.  Our  review  here  will  be 
necessarily  very  incomplete  but  will  touch  on  those  aspects 
of  particular  relevance  to  magnetic  bubble  applications. 
These  include  the  relationships  of  FMR  frequency  and  applied 
field  to  anisotropy  and  of  the  FMR  linewidth  to  mobility  and 
other  dynamic  properties  of  domains. 

The  most  useful  experimental  configuration  for  the 
present  purposes  is  that  generally  referred  to  as  perpen¬ 
dicular  resonance.  In  it  a  static  magnetic  field  H  is 
applied  perpendicular  to  the  film  plane  and  an  oscillating 
field  (generally  in  the  microwave  range)  is  applied  in  the 
plane  of  the  film.  The  condition  for  the  lowest  frequency 
mode  in  this  case  is  given  by 

~  =  H  -4ttMs  +  Ha  (12-1) 
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(Ref . 162 , 163 ,164 ,  35,  and  13).  Here  w  is  the  resonance 
frequency  and  y,  the  gyromagnetic  ratio.  This  equation  can 
be  pictured  in  the  following  elementary  way:  The  resonance 
condition  for  an  ion  is  simply,  w  =  YHtota^  where 
made  up  of  three  terms:  1)  The  externally  applied  field,  H, 
assumed  large  enough  to  saturate  the  sample,  2)  the  demagnet¬ 
izing  field,  -  4nMg ,  which  is  the  internal  magnetic  field 
generated  by  the  free  magnetic  poles  on  the  surfaces  of  the 
films  and  3)  the  anisotropy  field,  Hft,  which,  by  its  very 
definition,  plays  the  role  of  an  additional  bias  field 
resisting  the  tendency  of  the  microwave  field  to  tip  the 
net  magnetization  away  from  the  film  normal.  The  full 
spectrum  of  higher  excitations  in  FMR  is  very  rich,  as  the 
above  references  show.  In  addition  there  are  a  number  of 
effects  of  sample  geometry  and  orientation  which  we  will 
not  discuss.  (See,  for  example  references  165andl66).  For 
the  films  of  interest  in  this  report  the  magnetization  is  low 
enough  that  these  higher  frequency  modes  are  not  strong. 

Heinz,  et  al  (Ref.  35)  have  pointed  out  that  equation 
(12-1)  is  directly  relevant  to  the  conditions  for  babble 
stability  and  mobility.  In  order  to  have  bubbles  with 
magnetization  normal  to  the  film  plane  the  condition 

H.  >  4ttM  (12-2) 

A  **  s 

must  be  satisfied  (Ref.  29).  This  can  be  read  directly  from 
the  FMR  results  as 

Ha  -  4ttM  =  £  -  H  >  0  (12-3) 

A  S  y 
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An  additional  consideration  is  that  increasing  anisotropy 
causes  a  reduced  domain  wall  width  and  with  it  decreased 
mobility  so  that  maintaining  2  £  Ha/4ttMs  £  8  is  desirable. 

FMR,  coupled  with  a  measurement  of  4ttMs  as  described 
earlier  and  published  values  of  y ,  provides  useful  know¬ 
ledge  of  H.  and  this  ratio.  Pearson,  et  al  (Ref.  95 ,  96 , 

n 

and  97)  have  reported  detailed  comparisons  of  anisotropy 
data  from  FMR  and  torque  measurements  for  the  garnets.  The 
phenomenon  of  growth  induced  anisotropy  in  the  garnets, 
which  has  given  the  crystal  grower  much  wider  latitude  in 
designing  new  bubble  materials,  was  investigated  using  FMR 
by  LeCraw,  et  al  (Ref. 167).  FMR  has  the  capability  of 
providing  4^  and  y  values  if  several  frequencies  or 
orientations  are  used  but  such  an  approach  is  quite  involved 
and  does  not  lend  itself  to  routine  characterization. 

It  is  worth  pointing  out  that  the  connection  between 
static  torque  magnetometer  measurements  and  FMR  is  quite 
close.  This  results  from  the  fact  that  FMR  detects  the 
condition  in  which  maximum  energy  can  be  pumped  into  the 
magnetic  system.  Since  such  a  system  is  associated  with  an 
angular  momentum  a  torque  is  required  to  bring  this  about. 
Thus  the  starting  vector  equation  for  the  FMR  calculation  is 

=  y  M  x  H  =  ( 12-4 ) 

dt 

where  x  is  the  torque.  As  demonstrated  by  Pearson  and 
Cooper  (Ref.  95  -  see  also  ref. 16 8)  for  Tb  doped  YIG,  the 
field  for  resonance  is  given  by 
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which  is  also  the  second  derivative  of  the  crystal  energy 
with  respect  to  angle.  Here  0  is  the  angle  which  the 
applied  field  makes  with  a  crystal  axis.  Unfortunately, 
such  a  simple  relationship  only  holds  for  spherical  samples, 
reducing  its  usefulness  in  the  present  context. 

The  connection  between  domain  wall  mobility  and  FMR 
linewidth  was  developed  largely  during  the  1950's  by  workers 
at  Bell  Telephone  Laboratories  (Ref .119 ,169,118,  andl20). 
Hagedorn  and  Gyorgy  (Ref. 122)  noticed  a  discrepancy  between 
damping  factors  derived  from  mobility  and  FMR  measurements 
for  high  mobility  materials,  and  developed  a  new  theory  for 
wall  damping  in  such  materials.  Their  prediction  of  velocity 
proportional  to  the  square  root  of  the  drive  field  has 
recently  been  verified  by  Calhoun,  et  al  (Ref. 115). 

Harper  and  Teale  (Ref. 125)  pursued  the  connection 
between  mobility  and  FMR  and,  by  making  some  simplifying 
assumptions,  have  deduced  the  equation 

aV2 

Mw  s  ~ah~  nrrm2  F  (ai'  °2)  (12-6) 

where  AH  is  the  linewidth;  to,  the  resonant  frequency;  t, 
the  relaxation  time  for  the  excited  states  of  the  crystal 
and  F  (a^ ,  )  a  function  of  crystalline  anisotropy  and  the 

orientations  of  Mg  on  both  sides  of  a  domain  wall.  It  is 
clear  that  a  number  of  measurements  would  be  necessary  to 
extract  numerical  values  of  mobility  from  such  an  approach. 
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So  far  as  we  are  aware,  no  one  is  attempting  to  do  that  on  a 
routine  basis.  However,  this  connection  has  been  of  great 
use  in  the  search  for  improved  magnetic  bubble  materials. 
Because  of  the  large  amount  of  work  done  by  Dillon  and 
others  (Ref . 170 ,171 ,172 ,173 ,  and  94)  linewidth  information 
on  many  of  the  rare-earth  garnets  and  rare-earth  doped 
YIG  are  readily  available.  This  information  provides 
valuable  guidance  in  choosing  rare-earths  necessary  to 
adjust  the  magnetic  and  lattice  parameters.  It  indicates, 
among  other  things,  that  the  large  linewidth  ions  Tb ,  Dy, 
and  Ho  are  to  be  avoided  in  applications  demanding  high 
mobility. 

The  apparatus  and  techniques  for  FMR  are  well  established 
and,  at  least  for  room  temperature  work,  largely  available  on 
an  off-the-shelf  basis  today.  Kip  (Ref. 174),  Soohoo  (Ref. 166) 
and  Seiden  (Ref. 175)  have  discussed  apparatus  for  this 
measurement.  Fig.  20  illustrates  a  simple  system  as  shown 
by  Seiden.  In  it  the  power  reflected  by  the  sample  is 
detected  in  the  crystal  detector  and,  after  amplification, 
recorded  directly.  For  weak  signals  such  as  those  arising 
from  dilute  systems  of  magnetic  ions  more  sophisticated 
methods  such  as  those  used  in  electron  paramagnetic  resonance 
are  available  (Ref .176). 
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Figure  20.  Diagram  of  simple  FMR  apparatus  for  microwave  frequencies. 
After  Seiden,  ref.  175. 


13.  ANISOTROPY 

The  magnetic  anisotropy,  as  mentioned  previously  in  this 
report,  is  represented  by  ^he  anisotropy  field  required  to 
pull  the  magnetization  into  the  hard. direction,  ,  or  the 
related  uniaxial  anisotropy  constant,  K  =  HaM  /2 .  While 
this  uniaxial  anisotropy,  favoring  a  magnetization  direction 
normal  to  the  plane  of  the  bubble  film  or  platelet,  is 
essential  for  the  bubble  device,  it  is  an  over-simplification 
to  ignore  other  forms  of  anisotropy  which  may  be  present.  We 
therefore  present  a  very  brief  introduction  to  the  subject 
here.  More  extensive  discussion  can  be  found  in  a  number  of 
references  (see  for  example,  ref. 177  andl78) . 

The  orthoferrites  have  an  orthorhombic  crystal  structure 

and  exhibit  a  very  high  anisotropy,  M  lying  along  the  c 

s 

direction  at  room  temperature  for  most  of  the  single  rare- 
earth  compositions.  (The  lattice  parameters  obey  the 
inequality  a  <  b  <  c.)  The  exception  is  SmFe03  for  which  the 
a  direction  is  preferred  (Ref.  90,179,  andl80).  Such  a  di¬ 
rection  is  also  preferred  below  the  reorientation  temperature 
in  other  orthoferrites  (Ref. 181  and  section  16).  Away  from 
this  temperature  the  anisotropy  field  has  been  estimated  at 
approximately  10^  Oe  (Ref. 182).  This  gives  rise  to  a  large 
wall  energy  and  the  large  bubbles  which  have  been  the  major 
drawback  of  orthoferrites  for  devices.  By  growing  mixed 
rare-earth  orthoferrites  with  Sm  it  is  possible  to  come  close 
enough  to  the  reorientation  temperature  to  reduce  this 
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anisotropy  and  the  bubble  size  to  some  extent  but  only  at  the 
expense  of  temperature  stability. 

The  garnets  are  basically  cubic  crystals  which,  in  the 
absence  of  other  factors,  exhibit  typical  cubic  dependence  of 
energy  upon  magnetization  direction, 

E  =  Ki  (a12a22+a22a32+a32a12)  +  ...  (13-1) 

Here  is  the  direction  cosine  between  the  magnetization 
direction  and  the  ith.  cubic  axis.  Pearson  (Ref.  96)  has 
measured  K1  and,  where  possible,  K2  for  most  of  the  single 
rare-earth  garnets.  Generally  is  negative  leading  to 
easy  directions  for  the  magnetization  of  <111>  .  (See  also 
ref.  64).  Because  of  the  high  symmetry  exhibited  by  this 
expression  garnets  were  not  at  first  looked  upon  as  likely 
prospects  for  bubble  devices.  It  has  proved  possible, 
however,  to  generate  sufficient  uniaxial  anisotropy  to 
markedly  change  those  prospects.  Both  stress,  generated  by 
polishing  or  mismatch  with  the  underlying  substrate  (Ref.  104 
and  3  )  and  particular  growth  direction  and  constituents  - 
termed  growth  induced  anisotropy  (Ref. 183,  64,184,  andl85)- 
have  been  shown  to  yield  the  necessary  uniaxiality.  This 
uniaxial  anisotropy  can  be  described  as  a  series  of  terms 
in  the  energy  once  again: 

E  =  Ku]sin2e  +  Ku2sin40  +  ...  (13-2) 

Here  ©  is  the  angle  between  the  magnetization  and  the  easy 
axis.  The  total  anisotropy  energy  of  a  garnet  in  general 
then  will  be  given  by  the  sum  of  equations  (13-1)  and  (13-2). 
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For  device  purposes  the  most  important  anisotropy  is  between 
the  film  normal  direction  and  the  lowest  energy  direction 
approximately  in  the  film  plane  since  it  is  this  difference 
which  will  determine  the  bubble  stability.  This  difference 
is  here  referred  to  as  Ku* 

13 . 1  Non-Optical  Methods 

Two  of  the  standard  methods  for  the  measurement  of 
magnetic  anisotropy  have  already  been  discussed  in  some 
detail.  These  are  torque  magnetometry  (Section  9)  and 
ferromagnetic  resonance  (Section  12) . 

One  other  recently  developed  method  deserves  mention 
here.  Flanders  and  Doyle  (Ref. 186  andl87)  have  developed 
a  rotating  sample  magnetometer  and  applied  it  to  the 
measurement  of  anisotropy  in  GdIG  platelets.  The  voltages 
induced  in  pickup  coils  near  the  rotating  sample  are  detected 
using  a  lock-in  amplifier  synchronized  to  the  rotation 
frequency  or  one  of  its  harmonics.  Reference  186  discusses  a 
number  of  sample  configurations  to  be  used  depending  upon  the 
sample  anisotropy  and  the  property  desired.  It  indicates  that 
a  rotating  sample  system  may  have  advantages  over  vibrating 
sample  or  coil  systems  of  ruggedness  and  simplicity  while 
still  maintaining  comparable  sensitivity.  Such  a  system 
should  be  applicable  to  epitaxial  films  on  paramagnetic  sub¬ 
strates  but  with  the  same  general  limitations  discussed  for 
torque  magnetometry  in  Section  9. 
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13.2  Magneto-Optical  Methods 

Kurtzig  and  Hagedorn  (Ref. 188)  have  described  a  visual 
technique  in  which  a  variable  field  is  applied  in  the  plane 
of  the  sample.  The  Faraday  contrast  between  adjacent  domains 
decreases  as  the  field  is  increased  and  the  field  at  which 
the  contrast  is  undetectable  is  taken  as  H  .  The  author 

n 

claims  10%  accuracy  for  this  method.  They  report  results 
for  both  CVD  and  LPE  grown  films  including  the  observation 
that  annealing  at  temperature  in  excess  of  1200 °C  removes 
essentially  all  of  this  non-cubic  anisotropy  of  LPE  films 
but  does  not  affect  the  (strain  induced)  non-cubic  aniso¬ 
tropy  of  CVD  films.  This  visual  method  is  not  entirely 
satisfactory  in  that  the  rotation  of  magnetization  into  the 
hard  direction  is  a  gradual  function  of  applied  field  with¬ 
out  a  well  defined  visual  cutoff.  m  addition  saturation  of 
the  sample  by  a  small  normal  component  of  the  applied  field 
also  leads  to  loss  of  contrast,  and  could  be  misinterpreted 
as  Ha,  although  it  is  true  that  in  this  case  the  two  domain 
types  will  not  maintain  equal  areas.  A  more  quantitative 
approach  is  desirable. 

Such  a  method  has  been  developed  by  Shumate  (Ref.  47 
and  40) .  (See  also  reference  39  ) .  it  involves  applying 
a  field  Hx  in  the  plane  of  the  sample,  once  again,  with  the 
addition  of  an  alternating  field  Hz  normal  to  the  sample. 

As  long  as  domain  walls  are  present  they  move  in  synchronism 
with  this  alternating  component  (providing  a  domain  wall 
motion  term  to  the  susceptibility)  and  their  motion  is 


111 


readily  detected  using  a  silicon  or  photomultiplier  detector 

and  lock-in  amplifier.  Once  the  walls  have  been  eliminated 

a  significant  spin  susceptibility  remains  (the  spins  being 

readily  turned  from  the  hard  direction)  which  gradually  falls 

to  zero  as  the  in-plane  field  is  increased  further.  It  is 

convenient  to  add  provision  for  a  DC  component  to  H  -  H  (DC) 

z  z 

(normal  to  the  sample)  by  use  of  a  suitable  power  supply 

for  the  coil  which  provides  the  alternating  field.  Such  a 

system  is  shown  schematically  in  Fig. 21a. 

At  this  point  there  are  at  least  two  ways  to  determine 

from  the  system  output.  In  the  approach  favored  by 

Shumate,  Hz (DC)  is  adjusted  to  yield  saturation  of  the 

sample.  This  eliminates  the  domain  walls  and  leads  to  zero 

signal  until  Hx  reaches  approximately  HA  when  the  spin 

susceptibility  reappears  as  shown  in  Figure  21,  curve 

The  peak  in  this  curve  is  then  identified  as  H  -  4ttM  . 

A  s ' 

when  consideration  is  given  to  the  shape  anisotropy  or 

demagnetizing  factor  of  these  thin  films.  The  field  at 

which  this  peak  occurs  is,  in  fact,  a  function  of  the 

strength  of  Hz(DC)  as  shown  in  the  figure.  Shumate  (Ref. 189) 

has  calculated  a  correction  curve  which  takes  account  of 

this  effect  and  should  be  available  shortly.  The  present 

authors,  as  an  alternative  to  this  correction,  adjust 

Hz (DC)  until  the  wall  susceptibility  contribution  extends 

to  the  maximum  H^.  This  approach  amounts  to  using  H  (DC) 

z 


to  maxe  small  adjustments  in  the  total  field  direction 
until  that  direction  coincides  with  the  crystalline  hard 
direction.  The  maximum  in-plane  field  which  still  yields 
wall  motion  is  then  taken  as  -  4^.  Figure  21  shows 
several  curves  taken  at  different  values  of  H  (DC)  by 
this  technique.  Three  or  more  such  curves  are  required 
to  pinpoint  Hfl  -  4r«s,  but  these  can  be  obtained  rapidly 
due  to  the  substantial  signal  available.  Partly  because 
of  uncertainties  in  interpretation  this  method  is  estimated 
to  yield  Hfl  to  +  10*  but  is  felt  to  be  the  best  routine 
method  for  its  determination  in  epitaxial  thin  films  at  the 
present  time.  If  the  work  of  Shumate  (Ref.  189)  removes  this 
interpretational  uncertainty,  accuracy  of  4  3*  seems 
quite  possible.  Shumate  (Ref.  40)  has  also  used  this  technique 
to  study  the  in-plane  anisotropy  of  magnetic  bubble  films, 
thereby  further  extending  the  usefulness  of  this  method. 

With  possible  minor  modifications  in  the  future  it  appears 
to  be  by  far  the  most  convenient  approach  to  the  study  of 
anisotropy  in  films  for  bubble  devices. 
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14.  COERCIVITY 


The  coercivity,  Hc,  is  the  field  in  excess  of  the 
equilibrium  value  which  is  required  to  cause  a  domain  wall 
to  move,  i.e.,  one  half  of  the  hysteresis  loop  (M  vs.  H) 
width.  As  such  it  represents  an  energy  dissipation  mechanism 
to  be  held  to  a  practical  minimum.  It  is  well  known 
that  in  most  magnetic  materials  the  hysteresis  loops  for 
small  field  excursions  are  not  ideal  parallelograms  in 
which  no  change  in  M  occurs  until  the  major  hysteresis  loop 
is  intersected.  Rather,  the  minor  loops  continue  to  exhibit 
decreasing  slope  near  H=0  as  the  extent  of  the  excursions 
is  reduced.  The  magnetic  bubble  materials  are  rather  well 
behaved  in  this  respect  but  such  effects  do  make  the  coercivity 
a  somewhat  nebulous  parameter. 

14 . 1  Full  Hysteresis  Loop  Methods 

Historically  the  coercivity  has  been  measured  as  the 
half  width  of  the  full  hysteresis  curve,  i.e.,  the  M  vs.  H 
loop  which  extends  to  saturation  in  both  positive  and  negative 
applied  field  directions.  Such  a  hysteresis  curve  has 
been  presented  in  Figure  16  for  a  bubble  garnet  and  similar 
curves  have  been  published  by  C.  D.  Mee  (Ref . 104andl05)  and 
others  (Ref . 190 , 191 ,  andl92) .  All  of  these  references 
employed  the  Faraday  effect  as  discussed  in  Section  7  and 
elsewhere  (Ref . 193 ,194 ,  andl03) .  In  addition,  several  of 
the  methods  discussed  for  the  measurement  of  magnetization 
(Section  9)  will  yield  hysteresis  curves  and  a  measure  of 


the  coercivity.  These  include  mutual  induction  effects 
(Ref. 195)/  torque  magnetometry  (Ref.  93)  and  vibrating 
sample  magnetometry  (Ref. 179).  The  ballistic  measurement 
of  the  applied  field  required  to  yield  zero  net  magnetization 
has  also  been  used  (Ref. 196).  However,  the  Faraday  effect 
remains  the  best  method  for  the  measurement  of  the  hysteresis 
loop  and,  from  it,  the  coercivity  for  thin,  optically 
transparent  samples . 


14.2  Microscopic  Coercivity  Methods 

Bobeck,  et  al  (Ref. 153)  have  proposed  a  technique 
which  is,  once  again,  closely  related  to  the  bubble  domain 
phenomenology.  They  make  use  of  the  repulsive  interaction 
between  bubbles  to  yield  a  measure  of  the  smallest  force 
required  to  move  a  bubble.  This  is  done  by  bringing  two 
bubbles  together  (by  means  of  a  magnetic  probe,  for  example) 
and  then  allowing  them  to  repel  to  the  relaxed  distance, 

1^2'  By  approximating  the  field  of  a  bubble  as  that  of  a 
dipole  at  the  interaction  distance  a  coercivity  is  derived  as 


H  3irr  3h 

c  o 

4ttM  82,../ 

s  12 


(14-1) 


where  rQ  is  the  bubble  radius.  The  difficulties  in  using 
this  technique  for  accurate  measurements  in  materials  where 
the  bubble  diameter  is  small  are  immediately  evident  from 
the  high  powers  to  which  the  microscopic  distances  are  raised. 
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Ignoring  the  relatively  small  uncertainty  in  h,  standard 
error  propagation  yields 


AH  /H 
c  c 


9  (Ar  /r  ) 
o'  o 


16  (AA12/A12) 


1/2 


(14-2) 


Thus  these  uncertainties  are  amplified  by  these  large 
exponents  making  the  technique  impractical  when  r  and  &12 
cannot  be  measured  to  high  precision. 


Another  approach  to  the  coercivity  on  a  microscopic 
scale  was  developed  by  Kurtzig  (Ref.  50)  and  discussed  in 
connection  with  defect  detection  in  Section  7.  It  is  based 
upon  the  generation  of  a  straight  wall  by  a  gradient  field 
and,  in  the  garnets,  suffers  from  the  difficulties  already 
described. 

The  Bobeck  reference  (Ref. 153)  mentions  another  approach 
to  coercivity  on  a  microscopic  scale  which  is  a  more  workable 
alternative  for  the  garnets.  If  any  of  the  bubble  translation 
methods  for  the  measurement  of  mobility  are  used  at  low 
pulsed  gradient  fields  a  limit  can  be  found  below  which 
no  bubble  motion  takes  place.  Thus,  for  example,  in 
Equation  11-12  the  coercivity  is  determined  by  extrapolating 
the  velocity  vs.  pulse  current  curve  to  v=0  (at  I  )  and 
then  solving  the  resulting  equation  to  yield 


H 
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(14-3) 


Such  a  method  is  far  less  sensitive  to  the  uncertainties  in 
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bubble  diameter  and  is  accomplished  at  the  same  time  that 
the  mobility  is  being  determined.  This  technique  for 
mobility  measurement  (Ref. 158),  when  used  with  the  superimposed 
ramp  field  in  order  to  maintain  the  field  at  the  bubble 
constant,  is  judged  to  be  too  complex  to  be  used  in  routine 
characterization,  unfortunately.  For  such  purposes  we 
recommend  the  following  method  for  the  determination  of 
coercivity  with  the  recommended  mobility  technique  as  a 
separate  measurement. 

14.3  Partial  Hysteresis  Loop  Method 

The  hysteresis  methods  discussed  above  are  somewhat  mis¬ 
leading  for  the  purposes  of  magnetic  bubble  devices  in  that 
the  loops,  when  carried  to  saturation,  are  influenced  by 
nucleation  processes  at  the  edge  of  the  film  and  by  local 
defects  which  may  not  be  in  the  area  of  interest.  Minor 
hysteresis  loops  which  do  not  progress  to  saturation  are 
a  better  measure  of  the  local  coercivity  but  the  low 
signal  to  noise  ratio  becomes  a  serious  limitation  on 
accuracy.  Luff  (Ref.  53)  Shaw,  et  al  (Ref.  39)  and 
Shumate  (Ref.  40)  have  used  an  AC  technique  which  takes 
advantage  of  the  narrow  band  capabilities  of  the  lock-m 
amplifier.  In  it  the  domain  walls  are  driven  by  an  alternating 
bias  field  and  the  modulated  Faraday  signal  detected  using 
a  photomultiplier  and  the  lock-in  amplifier.  A  schematic 
diagram  of  the  apparatus  for  this  measurement  is  shown 
in  Figure  22.  The  resulting  signal  is  plotted  versus  the 
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Lock-in  Amplifier  Output 


Figure  22a 


Figure  22b 


Figure  22.  Coercivity  system 
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alternating  field  amplitude.  If  the  hysteresis  loops  were 

ideal  parallelograms  in  which  no  change  in  M  occurred  across 

the  top  and  bottom  of  the  loop,  then  the  lock-in  amplifier 

signal  would  fall  to  zero  at  an  alternating  field  amplitude 

of  Hc  and  remain  zero  for  all  smaller  amplitudes.  The 

typical  data  of  Figure  22  show  that  this  is  not  quite 

fulfilled  but  that  a  linear  portion  of  the  curve  is  apparent 

which  can  be  extrapolated  to  zero  lock-in  signal  to  yield 

H  .  The  short  vertical  lines  on  this  figure  are  the  data 
c 

points  taken  by  setting  the  recorder  pen  down  for  approximately 
10  time  constants  of  the  lock-in  amplifier  filter  circuit. 

The  "no  sample"  curve  is  similar  data  taken  with  no  sample 
in  the  microscope  and  results  from  Faraday  rotation  in  the 
microscope  optics. (Ref.  39).  Such  an  effect  can  be  avoided 
by  using  laser  illumination  and  avoiding  the  use  of  a 
microscope  entirely  or  by  using  long  working  distance 
optics  and  a  small  coil  (Ref. 155),  but  cross  checks  have 
indicated  the  present  technique  is  as  accurate  and  it  is 
more  convenient  to  carry  it  through  while  the  sample  is  in 
the  microscope  for  other  measurements. 
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15 .  MAGNETOSTRICTION  COEFFICIENTS 

In  single  crystal  magnetic  materials  the  magnetocrystalline 
anisotropy  energy  depends  on  the  state  of  strain  of  the  crystal 
lattice.  When  the  lattice  deforms  under  the  influence  of  a 
given  stress,  the  increase  in  elastic  energy  induces  a  change 
in  the  anisotropy  energy  so  as  to  minimize  the  total  energy  of 
the  system.  The  material  parameters  which  characterize  this 
change  in  anisotropy  energy  are  the  magnetostriction  coefficients, 

^100  anC*  *111*  T^ie  ma9netostr^ct^on  phenomenon  is  reciprocal 
in  the  sense  that  the  dimensions  of  the  crystal  depend  on  the 
direction  of  magnetization  relative  to  the  crystal  axes.  For 
materials  with  cubic  symmetry,  such  as  the  ferrimagnetic  garnets, 
the  magnetically  induced  dimensional  change  can  be  expressed  as 
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SL 
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+  3  *111  <0lVie2  +  a2a3e2B3  +  VlW'  <15-1) 


where  a^,  a^,  are  the  direction  cosines  of  the  magnetization 

direction  referred  to  the  cubic  axes,  and  0^,  02»  03  are  the 

A£ 

direction  cosines  of  the  direction  xn  which  -j  is  measured. 

The  change  in  elastic  energy  resulting  from  the  crystal 
dimensional  change  — j-  may  be  regarded  as  a  first  order  correction 
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to  the  anisotropy  energy,  and  as  such,  adds  a  term  to 
Equation  13-1.  Physically,  the  significance  of  Equation  15-1 
may  be  considered  in  the  following  light:  It  may  be 
energetically  favorable  for  the  crystal  to  deform  in  a  given 
direction,  if,  as  a  result,  the  anisotropy  energy  is  lowered  by 
more  than  the  elastic  energy  is  increased. 

Magnetostriction  in  magnetic  bubble  materials  is  important 
for  two  primary  reasons:  1)  In  thin  layers  of  ferrimagnetic 
garnets,  magnetostriction  has  been  utilized  extensively  to 
induce  sufficient  uniaxial  anisotropy  for  the  material  to 
support  magnetic  bubble  domains  (References  3,  35,  197,  and  198); 
2)  local  variations  in  the  state  of  strain  of  the  crystal 
lattice  (for  example,  near  a  dislocation,  lead  to  local 
inhomogeneities  in  magnetocrystalline  anisotropy,  and  inter¬ 
fere  with  bubble  propagation. 

The  earliest  technique  used  to  measure  the  magnetostriction 
coefficients  of  ferrimagnetic  garnets  was  employed  by  Callen 
et  al  (Referencel99)  to  measure  the  constants  of  single  crystal 
YIG.  The  method  utilizes  a  strain  gauge  apparatus  first 
employed  by  Goldman  (Reference200)  to  measure  magnetostriction 
constants  of  polycrystalline  rods  and  strips.  In  this  technique 
a  thin  metallic  foil  with  attached  leads  is  bonded  to  a  magnetic 


sample.  A  similar  device  is  bonded  to  a  dummy  ceramic  block 
in  close  proximity  to  the  first  so  that  both  will  experience 
the  same  temperature  and  magnetic  field  environment  during 
measurement.  The  resistance  of  the  foil  gauges  as  a  function 
of  strain  has  previously  been  calibrated  (usually  by  the 
manufacturer,  e.g.,  Budd  Company,  who  then  quotes  a  gauge 
factor)  so  that,  when  the  active  and  the  dummy  gauges  are 
connected  as  two  legs  of  a  D.  C.  Wheatstone  bridge,  resistance 

-7 

changes  due  to  sample  strains  as  small  as  10  can  be  detected. 
This  strain  gauge  technique  was  also  used  by  Iida  (Reference20l) 
to  measure  the  magnetostriction  coefficients  of  Tb3Fe5012  and 
Eu  Fe  0  ,  by  Clark  (Reference202)  to  measure  the  constants 

of  dysprosium,  holmium  and  erbium  iron  garnets,  by  Iida 
(Reference203 )  to  measure  the  constants  of  all  of  the  rare 
earth  iron  garnets,  by  Belov,  et  al  (Reference204)  to  measure 
the  constants  of  thulium  orthoferrite  in  the  reorientation 
temperature  range,  and  by  Gyorgy,  et  al  (Reference205)  to 
measure  several  partially  substituted  Sm  orthoferrites.  The 
data  of  Clark  and  Iida  are  not  in  precise  agreement  at  room 
temperature,  the  absolute  magnitude  of  Clark's  magnetostriction 
values  being  higher  in  general.  Current  magnetic  bubble 
literature  (Reference  64)  quotes  Iida's  data  almost  exclusively, 
but  this  may  be  due  more  to  a  consideration  of  where  the  work 
was  performed  (Bell  Laboratories)  than  to  accuracy  considerations 
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Callen,  et  al  (Reference  199)  points  out  that  magneto¬ 
resistance  effects  in  the  metallic  foil  at  low  temperatures 
makes  strain  gauge  measurements  taken  at  liquid  nitrogen 
temperature  somewhat  in  error.  These  authors  proposed  an 
alternate  measurement  technique  in  which  the  sample  was 
sandwiched  between  a  metal  support  plate  and  an  aluminum 
plated  silica  disk  which  served  as  one  plate  of  a  capacitor. 

A  Be-Cu  ring  placed  around  the  sample  served  as  the  other 
capacitor  plate.  The  whole  capacitor  assembly  was  made  part  of 
a  free-running  Hartley  oscillator  whose  frequency  was  monitored. 
For  small  changes  in  sample  length  induced  by  rotating  the 
assembly  in  a  magnetic  field,  the  magnetostriction  coefficients 

are  proportional  to  frequency  shift  of  the  oscillator.  A 
-7 

strain  of  10  produced  a  frequency  shift  of  about  10  Hz  when 
operated  at  5  MHz.  This  method  eliminates  magnetoresistance 
effects  as  well  as  effects  due  to  changes  in  gauge  factor  as  a 
function  of  temperature,  but  the  distributed  capacity  of  the 
oscillator  must  be  known  accurately,  and  there  are  the  additional 
bonding  problems  which  arise  because  of  the  different  thermal 
expansion  coefficients  of  the  sample,  support,  and  capacitor 
plates. 
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Another  technique  for  measuring  the  magnetostriction 
coefficients,  in  a  sense  the  reciprocal  of  the  methods 
discussed  thus  far,  was  employed  by  Petrakovskii  (Reference206 ) 
to  measure  the  constants  of  Ga-YIG.  This  method  consists  of 
mounting  a  single  crystal  spherical  sample  in  a  press  which 
was  integral  with  a  specially  designed  waveguide  section. 

A  known  amount  of  stress  could  be  applied  to  the  sample  by 
loading  the  press  with  weights.  Small  shifts  in  the  ferro¬ 
magnetic  resonance  peak  which  resulted  from  the  uniaxial 
compressive  stress  were  detected  by  an  X-band  spectrometer 
(adapted  to  fit  this  purpose)  which  could  resolve  shifts  of  as 
little  as  0.005%  in  the  location  of  the  peak.  The  advantages 
of  this  technique  are  that  only  relatively  small  samples  are 
required  («  1  mm) ,  the  method  is  relatively  simple  to 
instrument,  the  measurement  results  in  accuracy  comparable  to 
the  strain  gauge  method,  and  the  technique  can  be  automated. 

On  the  other  hand,  only  those  materials  can  be  measured  which 
exhibit  a  narrow  ferromagnetic  resonance  line  (which  many  of 
the  ferrimagnetic  bubble  garnets  do  not) . 

All  of  the  techniques  discussed  are  applicable  only  to 
bulk  single  crystal  materials,  while  most  present  day  bubble 
materials  are  thin  epitaxially  grown  films  which  are  being 
stressed  by  the  substrate  on  which  they  are  grown.  The  stress 
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level  is  generally  too  poorly  known  to  obtain  accurate  magneto¬ 
striction  constant  measurements  from  the  shift  in  the  ferro¬ 
magnetic  resonance  line  as  a  function  of  stress.  Added  to 
this  problem  is  the  fact  that  most  useful  garnet  materials 
contain  relatively  complex  mixtures  of  rare  earth  and  yttrium 
ions  as  well  as  gallium  ions  partially  substituted  for  iron  on 
tetrahedral  sublattice  sites.  These  compositions  would  have 
to  be  grown  as  bulk  single  crystals  to  enable  measurement  of 
magnetostriction  constants  by  the  strain  gauge  technique. 

Giess,  et  al  (Reference 19 8)  have  discussed  a  few  of  the 
problems  involved  in  attempting  to  use  the  bulk  single  crystal 
iron  garnet  magnetostriction  coefficients  to  interpret 
observations  of  thin  film  gallium  substituted  rare  earth  iron 
garnets.  As  these  authors  point  out,  what  is  normally  done  is 
to  assume  the  magnetostriction  coefficients  in  a  gallium 
substituted  iron  garnet  are  not  diminished  by  gallium,  and  that 
for  a  mixed  rare  earth  garnet  the  magnetostriction  coefficients 
are  a  linear  combination  of  those  of  the  single  ion  rare  earth 
garnets.  The  latter  assumption  probably  isn't  bad,  but  the 
former  is  almost  certainly  in  error.  There  really  is  no 
accurate  way  to  measure  magnetostriction  constants  in  thin 
films  of  bubble  garnets.  Measurements  using  FMR  or  the  strain 
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gauge  technique  should  be  performed  on  bulk  single  crystals 
having  compositions  identical  to  those  of  the  bubble  garnet 
films  for  which  the  magnetostriction  coefficients  are  desired. 
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16.  TEMPERATURE  EFFECTS 

The  temperature  characteristics  of  bubble  materials 

are  relevant  for  a  variety  of  reasons,  the  most  important 

of  which  relates  to  the  eventual  device  application  in  „hich 

they  may  be  used.  The  bubble  diameter  is  a  function  of  all 

the  magnetic  parameters  which  characterise  a  sample  of 

bubble  material,  and  all  of  these  parameters  are  functions 

of  temperature.  if  the  operating  temperature  of  a  bubble 

device  varies,  then,  unless  careful  precautions  have  been 

taken,  so  will  the  bubble  diameter.  I„  the  worst  case 

the  change  in  bubble  diameter  with  temperature  will  result 

rn  bubble  collapse  or  runout  at  some  point  in  the  bubble 
device. 

16.1  Reorientation  Temperature  Range 

In  the  rare  earth  orthoferrites,  the  easy  axis  of  magneti- 
ration  undergoes  a  continuous  reorientation  from  the  crystal- 
g  aphic  c  to  the  a  axis  in  a  certain  region  of  temperature 
known  as  the  reorientation  temperature  range.  This  fact  along 
with  the  non-alterable  magnetisation  and  high  anisotropy 
of  the  orthoferrites  eliminated  them  from  consideration  as 
useful  bubble  domain  materials.  Because  of  the  drastic 
nature  of  the  change  in  anisotropy  with  temperature  in 
the  reorientation  range,  several  methods  could  be  used 
to  measure  its  onset  and  its  temperature  extent.  Sherwood, 
et  al  (Ref. 182)  used  a  torque  magnetometer  in  a  variable 
temperature  mode  to  detect  the  reorientation  temperature 
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range.  Gyorgv ,  et  al  (Ref. 207)  used  both  variable  temperature 
torque  and  vibrating  sample  magnetometers  for  reorientation 
range  measurements  in  Sm  and  Tm  orthoferrites.  Rossol  (Ref. 208) 
made  an  exhaustive  study  of  the  reorientation  temperature  range 
in  thulium  orthoferrite  by  making  visual  measurements  of 
domain  size  and  then  using  the  theory  of  Malek  and 
Kambersky  (Ref.  61)  to  obtain  wall  energy  as  a  function  of 
temperature.  All  of  these  measurement  methods  have  previously 
been  described  so  no  further  discussion  will  be  given  here. 

16.2  Neel  and  Compensation  Temperatures 

The  temperature  at  which  the  transition  from  ferro¬ 
magnetic  to  paramagnetic  behavior  occurs  in  magnetic  materials 
is  known  as  the  Curie  temperature.  In  ferrimagnetic 
materials  this  temperature  is  >’so  known  as  the  Neel 
temperature.  The  rare  earth  garnets  exhibit  an  additional 
significant  temperature  dependent  effect.  The  rare  earth 
and  octahedral  iron  sublattices  align  antiparallel  to  the 
tetrahedral  iron  sublattice,  and  the  various  sublattices 
have  different  temperature  dependences.  As  a  result,  most 
of  the  garnets  experience  a  temperature  where  the  oppositely 
directed  sublattice  magnetizations  cancel,  resulting  in 
zero  net  moment.  This  temperature  is  known  as  the  compensation 

temperature . 

The  Neel  temperature  of  orthoferrites  and  garnets  has 
been  measured  by  a  variety  of  techniques  nearly  all  of  which 
involve  monitoring  the  magnetization  as  the  temperature  is 
increased  through  the  Neel  temperature.  The  most  cumbersome 
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in  terms  of  experimental  apparatus  is  the  neutron  diffraction 
technique  used  by  Koehler,  et  al  (Ref. 209)  to  measure  the 
Neel  temperature  of  polycrystalline  Nd,  Ho  and  Er  orthoferrites 
These  authors  discuss  the  design  and  operation  of  their 
apparatus  in  an  earlier  article  (Ref. 210).  Treves  (Ref.  92) 
has  done  Neel  temperature  measurements  on  most  of  the 
orthoferrites  by  torque  magnetometer  and  Mossbauer  spectroscopy 
Eibschiitz ,  et  al  (Ref. 211)  have  also  measured  Neei  temperatures 
in  all  the  orthoferrites  by  means  of  the  Mossbauer  technique 
with  an  accuracy  estimated  to  be  +  1°K.  Eibschiitz,  et  al 
(Ref. 212)  had  previously  measured  the  Neel  temperatures  of 
the  orthoferrites  by  differential  thermal  analysis  which 
relies  on  a  specific  heat  anomaly  at  the  Neel  temperature. 

Using  NiO  as  a  reference,  they  claimed  an  accuracy  of 
+  2%  for  those  Neel  temperature  determinations. 

The  above  references  represent  some  of  the  more 
unusual  approaches  to  the  measurement  of  Neel  and  compensation 
temperatures.  In  fact,  most  measurements  of  these  properties 
in  f errimagnetic  garnets  have  used  torque  magnetometer 
and  ferrimagnetic  resonance  measurements  as  functions  of 
temperature  to  determine  Neel  and  compensation  temperatures 
in  bulk  single  or  polycrystalline  garnet  materials.  The 
reader  is  referred  to  the  excellent  bibliography  following 
chapter  3  in  von  Aulock's  Handbook  of  Microwave  Ferrite 
Materials  ( Re f . 2 1 3 )  . 

The  technique  used  at  Monsanto  to  measure  the  Neel 
temperature  is  a  magneto-optical  one  similar  to  that 
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described  by  Kurtzig,  et  al  (Ref.  65).  This  method  involves 
the  modulation  of  the  Faraday  effect  by  oscillating  domain 
walls  which  are  being  driven  by  an  AC  bias  field.  The 
garnet  sample  is  situated  in  a  heating  stage  and  as  the 
stage  temperature  is  raised  to  the  Neel  temperature  the 
Faraday  susceptibility  drops  sharply  to  zero,  allowing  a 
measurement  precision  of  +  1°K.  This  method  is  rapid,  simple 
and  accurate  and  should  become  the  most  widely  used  technique 
for  determining  Neel  temperature  in  the  garnets. 

A  somewhat  similar  but  less  accurate  technique  is 
to  observe  the  domain  pattern  visually  as  the  sample 
temperature  is  raised  to  the  Neel  temperature.  In  the 
vicinity  of  the  Neel  temperature  the  domains  will  disappear. 
Unfortunately,  the  contrast  is  low  near  the  Neel  temperature 
and  the  exact  temperature  of  the  ferrimagnetic  to  para¬ 
magnetic  transition  is  difficult  to  determine  in  this  way. 
Accuracies  of  no  more  than  +  5%  should  be  expected  with 

this  method. 

Epitaxial  layers  with  a  highly  defected  film-substrate 
interface  will  spontaneously  demagnetize  in  an  array  of 
bubbles  when  brought  down  from  above  the  Neel  temperature 
as  domains  nucleate  at  every  defect  as  well  as  at  the  sample 
edge.  This  fact  has  been  used  by  Hsu  (Ref. 214)  to  find  the 
Neel  temperature  in  F.uEr  iron  gallium  garnet  films.  A 
film  previously  saturated  at  room  temperature  is  heated 
to  near  the  Neel  termperature ,  then  cooled  and  the  type 


of  domain  pattern  (bubbles  or  strips)  is  noted.  This 
procedure  is  repeated  in  small  temperature  increments  until 
the  Neel  temperature  is  reached.  Unfortunately  not  only 
is  this  method  tedious,  but  in  high  quality  samples  of  the 
type  required  in  bubble  devices,  the  film-substrate 
interface  is  not  defected,  and  the  layer  does  not  spontaneously 
demagnetize  into  an  array  of  bubbles.  Therefore  this 
technique  is  not  recommended. 

16.3  Temperature  Variation 

It  is  important  to  know  how  the  magnetization  and/or 
bubble  diameter  varies  with  device  operating  temperature, 
as  pointed  out  earlier  in  this  section.  The  functional 
dependence  of  bubble  diameter  on  temperature  has  been 
observed  visually  in  most  cases  (Ref . 208and215) .  Magnetization 
as  a  function  of  temperature  can  be  obtained  from  any  of 
the  techniques  of  section  9  applied  at  the  temperatures 
desired.  For  the  magneto-optical  techniques,  heating  stages 
for  microscopes  are  available  from  several  commercial 
sources  (a  Leitz  350°C  stage  is  used  at  Monsanto)  which  can 
be  useful  above  room  temperature  in  this  application. 

A  better  approach  over  somewhat  more  limited  range  is  to 
employ  a  thermoelectric  heating/cooling  device.  This 
device  has  several  advantages:  the  new  equilibrium 
temperature  is  attained  rapidly,  both  heating  and  cooling 
are  accomplished  electrically,  and  it  is  economical  to 
instrument.  The  chief  disadvantage  is  that  temperature 
excursions  are  limted  to  about  +  50 °C  in  the  region  of 
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room  temperature.  In  most  cases,  knowledge  of  bubble 
bias  field  and  diameter  variations  in  this  temperature 
range  along  with  Neel  and/or  compensation  point  data  is 
enough  to  evaluate  the  material  in  question  in  light  of 
its  eventual  performance  in  a  device  related  environment. 
Temperature  dependence  of  several  of  the  other  properties 
of  magnetic  bubble  films,  including  characteristic  length, 
domain  wall  energy,  coercivity,  and  mobility,  are  also 
measurable  with  such  a  system. 

16.4  Thermal  Conductivity 

One  additional  thermal  characteristic,  thermal 
conductivity,  is  worth  mentioning  briefly  since  it  may 
become  important  in  bubble  device  packaging.  Luthi  (Ref. 216) 
measured  the  thermal  conductivity  of  single  crystal  YIG 
at  liquid  helium  temperatures  using  a  standard  thermal 
conductivity  apparatus  which  consisted  of  carbon  resistors 
for  temperature  measurement  cemented  to  the  rectangular 
YIG  specimen  and  a  Dahlohm  resistor  cemented  to  the 
specimen  for  heating.  One  end  of  the  specimen  was  cemented 
to  a  copper  bar  which  was  immersed  in  liquid  helium. 

Luthi' s  results  indicate  a  thermal  conductivity  varying  from 
about  6mw/cm°K  at  1.5°K  up  to  about  50mw/cm°K  at  5°K.  An 
absolute  accuracy  of  20%  is  claimed  for  these  data.  Douglass 
(Ref. 217),  using  a  nearly  identical  apparatus,  measured  the 
thermal  conductivity  of  YIG  over  a  temperature  range  from 
0.4  to  20°K  in  the  presence  of  magnetic  fields  of  from  0  to 
20K  Oe.  His  data  indicates  a  maximum  in  thermal  conductivity 
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of  approximately  1000  mw/cm°K  in  the  vicinity  of  20 °K. 
Douglass'  data  agree  well  with  those  of  Luthi  (Ref. 216) 
in  the  common  range.  Slack  and  Oliver  (Ref. 218)  have 
recently  published  results  for  several  non-magnetic  A1 
and  Ga  garnets  as  well  as  YIG  for  temperatures  between 
2  K  and  300 °K.  These  authors  quote  a  room  temperature  thermal 
conductivity  for  YIG  of  0.09  W/cm°K  and  an  accuracy  of 
about  10%  using  a  refined  version  of  the  same  apparatus  dis¬ 
cussed  above  (Ref. 219).  In  all  likelihood,  this  parameter 
will  vary  little  from  one  ferrimagnetic  garnet  to  another, 
obviating  the  need  for  its  routine  measurement.  Sufficient 
accuracy  should  be  obtained  in  most  instances  if  the  value 
of  thermal  conductivity  of  YIG  is  used  in  calculations 
involving  other  garnets. 
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17.  OPTICAL  PROPERTIES 

Optical  methods  have  been  discussed  as  possible  means 
of  detection  or  "read-out"  in  magnetic  bubble  memories 
(Ref .220and221) .  As  Strauss  concludes,  optical  readout 
is  probably  not  the  most  attractive  method,  particularly 
with  the  advent  of  the  small  bubbles  available  in  garnets. 
Nevertheless,  for  the  person  involved  in  characterizing 
such  materials  optical  methods  are  extremely  convenient  and 
because  of  this,  these  properties  take  on  considerable 
interest.  In  particular,  the  index  of  refraction  enters  the 
thickness  measurement  and  the  Faraday  rotation  is  central 
to  all  of  the  optical  study  of  domain  dimensions  and 
behavior.  Other  optical  properties  ,  such  as  the  rare  earth 
spectra,  impurity  and  lattice  mode  spectra,  Kerr  effect 
and  Zeeman  and  Stark  spectroscopy,  possess  some  potential 
usefulness  and  a  great  deal  of  basic  literature,  but 
are  peripheral  to  the  central  theme  of  this  reDort.  As  such, 
they  will  be  given  only  brief  consideration.  A  third  class 
of  optical  properties  --  those  relating  to  the  instruments 
used  in  these  measurements  are  very  important  to  the  worker 
in  this  field  but  outside  of  the  scope  of  this  report. 

Some  references  of  a  review  nature  may  be  helpful,  however, 
and  are  presented  here:  Monochromators  -  (Ref .222),  micro¬ 
scopes  -  (Ref . 223) ,  light  sources  (Ref. 224),  fiber 
optics  (Ref . 225and226)  and  optical  detectors  (Ref. 227). 
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17. 1  Refractive  Index 

The  first  measurements  of  the  index  of  refraction  of 
the  garnets  involved  the  measurement  of  reflectivity 
(Ref. 228).  Solving  equation  5-2  for  n  yields 


n 


1  +  JR 

=  i  -  yir 


(17-1) 


where  R  is  the  fraction  of  the  incident  light  intensity 
which  is  reflected  and  the  equation  considers  only  one 
sample  surface.  In  this  reference  natural  single  crystal 
faces  were  used  and  later  surfaces  presumably  did  not 
contribute  to  the  reflected  beam.  The  resulting  index  for 
YIG  and  ErIG  was  quoted  at  2.2  +  0.2  in  the  0.6  to  l.Oym 
wavelength  range.  Geometric  and  surface  quality  effects 
tend  to  limit  the  accuracy  of  this  approach  and,  in  addition, 
there  is  always  the  suspicion  that  the  surface  layer  from 
which  reflection  takes  place  may  not  be  representative 
of  the  bulk  material.  For  the  thin  samples  of  current 
interest  here,  there  are  the  additional  difficulties  of 
reflection  and  interference  from  the  other  surfaces  present. 

Johnson  and  Walton  (Ref.  41)  measured  the  infrared 
index  for  YIG  and  Yx  ^Gd-j^  5Fe5C>12  by  forming  prisms  from 
flux-grown  C3 /stals  and  measuring  the  deviation  of  the  beam 
through  them.  As  discussed  in  optics  toxts  (Ref. 229)  one 
normally  adjusts  the  prism  for  minimum  deviation  of  the 
refracted  beam  in  which  case 


n  = 


Sin  (4min  +  “l /2 


Sin  a/2 


(17-2) 


where  6  .  is  the  minimum  angle  through  which  the  incoming 
min 

beam  is  bent  by  the  prism  and  a  is  the  angle  between  the 
two  faces  of  the  prism.  These  authors  present  values  of  n 
between  1.4pm  (n=2.209)  and  5.5pm  (n=2.103)  with  a  mean 
deviation  of  0.003.  No  detectable  difference  was  found 
between  the  two  garnet  compositions. 

Grant  (Ref.  42)  published  reflectivity  data  for  YIG 

O  O 

and  YGG  over  the  wavelength  range  6000A  to  1700A.  Due  to 
surface  and  alignment  difficulties,  however,  his  absolute 
values  are  uncertain  to  ~  30%  so  that  the  index  values 
resulting  from  Equation  17-1  are  not  directly  useful. 
Unfortunately,  the  wavelength  range  covered  by  Johnson 
and  Walton  does  not  permit  a  sufficiently  reliable  extrap¬ 
olation  into  Grant's  range  to  fix  his  absolute  values 
much  better  although  it  definitely  indicates  that  they 
are  too  low.  When  used  with  the  index  values  of  Roman, 
et  al  at  6328A  (see  below) ,  the  data  of  Grant  do  serve  to 
indicate  the  further  growth  in  n  at  shorter  wavelengths. 

Kahn,  et  al  (Ref. 2 30)  calculated  the  real  and  imaginary 
parts  of  the  index  of  refraction  using  Grant's  reflectivity 
results  and  their  own  Kerr  effect  values.  The  absolute 
values  are  subject  to  the  uncertainties  already  mentioned 
in  connection  with  Grant's  data. 

Grismore  and  Rhodes  (Ref. 231)  have  reported  reflectivity 
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measurements  on  Gd  and  Tb  iron  garnets  in  the  energy 
range  2  to  5.5eV  {6200  to  2260A  wavelength).  They  employed 
light  polarized  in  the  plane  of  incidence  and  made  measure¬ 
ments  at  two  angles  of  incidence,  <{>,  in  order  to  determine 
both  the  real  part,  n,  of  the  index  of  refraction  and  the 
imaginary  part,  k,  (sometimes  called  the  extinction 
coefficient)  which  becomes  significant  in  this  strongly 
absorbing  range.  For  this  polarization  the  reflectivity 
is  given  by 


r  =  (n  -  sec  <t>)2  +  k2 

(n  +  sec  <J>)  2  +  k2  <17-3> 

The  computed  values  of  n  for  TbIG  show  some  structure 
but  are  approximately  2.15  +  0.05  between  2  and  3eV. 
Unfortunately  they  did  not  compute  n  and  k  for  GdIG  below 
3eV  and,  indeed,  it  would  appear  that  the  very  low 
reflectivities  measured  at  <f>  =  70°  for  this  compound 
would  require  unreasonably  large  n  values  at  2eV  (~  2.7). 
Neglecting  absorption  and  using  only  their  value  of  R 
at  <p  =  20°,  we  compute  n  =  2.42  at  2eV.  Their  curve 
indicates  n  =  2.52  at  3eV  making  our  value  at  lower  energy 
appear  reasonable.  This  difference  in  n  between  TbIG  and 
GdIG  indicates  that  significant  errors  may  result  from 
assuming  a  constant  n  for  all  garnet  compositions.  A 
systematic  study  of  the  variations  with  composition  would 
be  of  considerable  value. 
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The  Brewster  angle  method  employed  by  Roman,  Spiwak, 
and  Baron  (Ref.  36)  makes  use  of  the  polarizing  effect  upon 
reflection  at  the  air-garnet  film  interface.  Theory  shows 
that  light  polarized  in  the  plane  of  incidence  has  zero 
reflectivity  at  a  certain  angle  of  incidence.  This 
angle,  ,  measured  from  the  normal  to  the  reflecting 
surface  is  known  as  Brewster's  angle  and  is  given  by  (Ref. 232) 

tan  4>b  =  n  (17-4) 

for  light  incident  from  vacuum  (or  air  to  ~  0.03%  accuracy) 
onto  a  surface  of  index  n.  Thus,  an  accurate  measurement 
of  the  angle  corresponding  to  this  zero  of  reflected  light 
intensity  yields  an  accurate  measure  of  n.  The  use  of  a 
laser,  with  its  highly  collimated  beam,  aids  in  improving 
the  accuracy  of  the  measurement.  For  the  case  of  films  on 
substrates  of  different  index  and  for  thin  samples  generally, 
the  measurement  is  complicated  by  the  presence  of  the 
reflected  beam  from  the  other  surfaces.  In  general  these 
will  add  a  component  which  will  not  go  to  zero  at  the 
Brewster  angle  and  may,  in  fact,  give  rise  to  interference 
effects.  These  will  go  to  zero,  of  course,  at  the  Brewster 
angle.  For  a  thick  film  in  the  visible 

(a  region  of  significant  absorption,  see  below) ,  the  proper 
angle  will  correspond  to  a  minimum  in  the  reflected 
intensity  of  the  correct  polarization. 

An  alternative  approach,  due  to  Abeles  (Ref . 233and234) 
takes  advantage  of  the  zero  reflectivity  and  the  fact  that 
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this  effect  is  reciprocal  (i.e.  also  occurs  for  light 
progressing  along  the  same  path  from  within  the  sample) . 

For  a  film  thin  enough  to  neglect  absorption,  the  angle  <f>D 
above  can  be  found  as  th  t  angle  at  which  the  substrate 
with  film  gives  the  same  reflected  intensity  as  a  similar 
bare  substrate.  Both  approaches  to  the  measurement  of  1 
the  Brewster  angle  are  subject  to  the  possible  criticism 
that  they  only  sample  the  surface  region.  However,  when 
used  in  combination  they  appear  to  yield  the  best  method 
for  the  measurement  of  refractive  index  of  magnetic 
bubble  materials. 

Current  refractive  index  measurements  at  Bell  Labs 
are  being  made  by  W.  J.  Tabor  on  bulk  samples  of  the  same 
compositions  as  those  of  interest  for  films  (Ref. 155). 

The  method  employed  is  very  closely  related  to  the  thickness 
measurements  of  films  already  discussed.  The  samples  are 
prepared  with  two  very  flat  opposite  faces  slightly  non¬ 
parallel  to  one  another.  The  interference  fringes  produced 
in  monochromatic  light  together  with  the  difference  in  thickness 
across  the  sample  are  sufficient  to  yield  the  index  of  refraction. 
For  N  interference  fringes  between  points  of  thickness 
difference,  Ah,  the  index  is  given  as 


n 


AN 

2Ah 


(17-5) 


This  method  is  also  applicable  to  films,  of  course,  when  Ah 
can  be  obtained  with  sufficient  accuracy  by  some  other  means. 
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Microscopic  measurements  on  broken  or  polished  edges  of  a 

tapered  thick  film  is  such  a  method.  For  typical  films, 

however,  the  Brewster  angle  method  appears  preferable. 

The  same  techniques  described  here  for  the  garnets  are 

applicable  to  the  orthoferrites  without  further  complication. 

Infrared  reflectivity  data  have  been  reported  by  Antonov, 

et  al  (Ref. 235)  who  analyzed  them  to  yield  n  =  2.35  +0.05 

at  X  =  2.0ym.  This  is  approximately  0.16  higher  than  YIG 

at  the  same  wavelength.  More  detailed  analysis  of  erbium 

orthoferrite  has  been  carried  out  using  the  Kerr  effect 

(see  below)  by  Jung,  (Ref. 236)  based  upon  theoretical  work 

of  Argyres  (Ref. 237)  and  others  (Ref. 238  and  239).  The 

results  are  presented  as  the  real  and  imaginary  parts, 

e  and  e~,  of  the  dielectric  constant  over  the  range  3600  to 
1  2 

7000A.  Analyzing  these  on  the  basis  of  the  Equations  (Ref. 240) 

?  2 

n-k  =  and  nk  =  e2 

o  ° 

one  finds  n  =  2.50  at  5000A  and  n  =  2.32  at  3900A. 
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17.2  Faraday  Effect 

In  the  case  of  the  orthoferrites  the  measurement  of 
Faraday  rotation  of  the  plane  of  polarization  is  complicated 
by  the  simultaneous  presence  of  birefringence  in  the  crystal. 
This  causes  the  plane  of  polarization  to  oscillate  as  the 
wave  progresses  through  the  crystal  rather  than  turning 
unidirectionally  as  it  would  under  the  influence  of  the 
Faraday  effect  alone.  The  analysis  of  this  situation  is 
sufficiently  complex  and  excellent  discussions  so  readily 
available  that  it  would  be  pointless  to  try  to  condense 
them  here.  Some  of  these  discussions  are  given  in 
references  241,242,243,244,  and245.  The  second  of  these 
indicates  that  all  of  the  rare  earth  orthoferrites,  with 
the  exception  of  SmFe03  (which  has  a  reorientation 
temperature  above  room  temperature)  follow  very  nearly  the 
same  curve  of  intrinsic  Faraday  rotation  vs.  wavelength 
in  the  range  0.55ym  <  X  <  1.8ym. 

The  garnets  are  basically  cubic  and  therefore  do 
not  exhibit  the  phenomenon  of  birefringence  in  the  usual 
sense.  However,  Dillon  (Ref . 246and247) ,  in  his  initial 
reports  of  optical  transparency  in  the  garnets,  noted  a 
magnetic  birefringence  for  light  traveling  perpendicular 

j 

to  the  direction  of  the  magnetization.  (See  also  ref. 248) 
While  this  effect  is  useful  for  pathologic  samples  with 
the  easy  axis  in  the  plane,  it  will  not  concern  us  for 
device  quality  magnetic  bubble  films.  Birefringence 


& 
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associated  with  growth  induced  anisotropy  has  also  recently 
been  reported  by  Dillon,  et  al  (Ref. 249). 

Dillon  (Ref. 246)  also  made  initial  measurements  of  the 
Faraday  rotation  of  several  iron  garnets  for  light 
traveling  along  the  magnetization  direction  and  reported 

O 

those  for  YIG  in  the  7700  to  5300A  range.  The  results 
show  some  broad  line  structure  on  a  generally  rising  curve 
as  A  decreases  (800  to  4000  degrees  /cm  over  this  range) . 

In  this  case  the  measurement  consists  simply  of  determining 
the  angle  between  the  extinction  condition  with  and  without 
the  sample  in  the  beam  (for  a  single  domain  area  of  the 
sample) .  A  separate  measurement  of  the  thickness  is 
required  to  yield  the  Faraday  rotatory  power  quoted  above. 
For  the  polished  platelets  used  by  Dillon  very  standard 
techniques  for  this  thickness  measurement  sufficed. 

Since  these  initial  reports  a  number  of  Faraday 
studies  in  the  garnets  have  been  reported.  Krinchik  and 
Chetkin  (Ref .228 and250)  have  reported  results  on  yttrium-, 
erbium-,  and  holmium  iron  garnets  at  77°K  and  300 °K  from 
less  than  lym  wavelength  to  ~  8ym.  These  show  a  wavelength 
independent  range  beyond  5ym  with  rising  rotation  and 
some  rare  earth  associated  structure  at  shorter  wavelengths. 
Johnson  and  Tebble  (Ref . 251and252)  have  reported  results 
in  the  1  to  5ym  wavelength  range  for  the  iron  garnets  of 
Y,  Y-Gd,  Gd,  Tm,  Dy  and  have  derived  g  factors  for  the 
ions  involved.  The  experimental  system  involved  recording 
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the  light  output  as  a  function  of  polarizer  angle  near 
extinction  as  described  by  Austin  (Ref. 253). 

LeCraw,  et  al  (Ref. 254)  were  concerned  with  the  ratio 
of  Faraday  rotation  to  light  absorption  in  the  near  infrared 
for  YIG,  determining  a  value  of  800  degrees/decibel  loss. 
MacDonald,  et  al  (Ref. 255)  indicate  that,  for  a  very  thin 
polycrystalline  GdIG  film,  this  value  is  less  than  unity. 
Dillon  (Ref. 256)  has  reviewed  the  properties  of  a  number 
of  magnetic  materials  for  potential  applications  as  Faraday 
rotators,  modulators,  lasers,  etc.  Almasi  (Ref. 221) 
has  considered  the  practicalities  of  magneto-optic  display 
and  bubble  sensing  systems  for  DyFeC>3 ,  YIG  and  Ga-YIG, 
and  PbFegAl^O^g.  Further  Faraday  effect  data  have  been 
given  for  Ga-YIG  by  Matthews,  et  al  (Ref. 257)  and  for 
GdIG  films  by  Sawatzky  and  Kay  (Ref. 258).  Sawatzky  and 
Horne  (Ref. 259)  have  described  a  highly  automated  system 
for  the  direct  recording  of  Faraday  rotation  and  coercivity 
which  employs  sampling  techniques  at  various  points  on 
the  optical  hysteresis  loop. 

17 . 3  Other  Optical  Properties 

The  Kerr  effect  -  the  transformation  of  plane  polarized 
light  upon  reflection  into  either  elliptically  polarized 
light  or  into  a  different  plane  of  polarization  -  is 
of  some  usefulness  in  observing  domain  patterns.  The  amount 
of  change  in  the  incident  light  is  very  much  smaller, 
however,  than  in  the  Faraday  effect  so  that  the  latter 
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will  doubtless  continue  to  be  used  wherever  possible.  A 
few  references  on  the  Kerr  effect  in  magnetic  bubble 
materials  are:  orthoferrites  -  References  260  and  230; 
garnets  -  References230and261.  These  papers  contain 
references  to  the  relevant  experimental  techniques. 

Optical  and  magneto-optical  absorption  studies  in  the 
orthoferrites  and  garnets  have  been  actively  pursued  since 
the  first  reports  of  optical  transparency.  They  provide 
considerable  insight  into  the  ionic  and  magnetic  structure 
of  these  crystals.  The  transmitted  light  intensity  is 
governed  by  the  law  I  =  IQ  exp (-ax)  where  IQ  is  the  incident 
intensity  (less  that  part  which  is  reflected) ,  I  is  the 
intensity  after  a  distance  x  is  traversed  in  the  sample, 
and  a  is  the  absorption  coefficient.  Except  for  line 
structure  associated  with  specific  rare  earth  ions,  all 
of  the  garnets  display  basically  the  same  wavelength 
dependence  as  shown  in  Figure  23.  The  apparatus  and 
methods  involved  are  generally  those  which  are  standard 
for  all  forms  of  spectrometry,  frequently  making  use  of 
commercially  available  spectrometers  with  facilities  for 
dual  beam  operation  and  chopped  light  beam  lock-in 
amplification  for  noise  suppression.  Some  representative 
publications  dealing  with  material  of  interest  to  the 
magnetic  bubble  worker  ere:  orthoferrites  -  References  262, 
263,  and264;  garnets  -  References  265,266,267,268,269,270, 
and  271. 
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Figure  23.  Optical  absorption  coef 
in  the  visible  and  near 
After  Wood,  ref.  270. 
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18.  ELECTRICAL  PROPERTIES 


The  magnetic  garnets  and  orthoferrites  are  highly  insula¬ 
ting  crystals  having  a  resistivity  in  excess  of  10"^  ohm-cm 
when  pure.  Wood  and  Remeika  (Ref. 272)  studied  the  effect  of 
Ca  doping  (leading  to  compensating  Fe4+  ions)  and  Si  doping 
(leading  to  Fe  ions)  on  optical  and  electrical  properties 
of  YIG.  In  both  cases  the  resistivity  dropped  to  ^  106  ohm 
cm  for  0.002  impurity  atoms  per  formula  and  leveled  off  at 

5 

%  10  ohm  cm  for  higher  concentrations.  In  addition,  they 
found  that  Fe4f  corresponds  to  p-type  conductivity  and  Fe2+ 
to  n-type  by  means  of  the  thermoelectric  effect.  Thus  the 
electrical  properties  of  garnets  are  demonstrated  to  be 
potentially  useful  as  a  semi-quantitative  check  of  purity. 
Similar  measurements  of  Sn  doping  of  Gd  and  Tm  orthoferrites 
(Ref. 26 4)  show  very  similar  pure  resistivities  with  decreases 

4 

to  'v.  10  ohm  cm  with  doping. 

Experimental  precautions  in  the  measurement  of  resis¬ 
tivity  of  such  insulating  materials  have  been  discussed  many 
times  in  the  literature  (Ref .273 ,274 ,  and275) .  They  consist 
primarily  of  1)  avoiding  shunting  paths  due  to  surface  effects 
and  instrumentation,  2)  avoiding  spurious  effects  of  contact 
resistance  and,  3)  avoiding  error  in  the  measurement  of  the 
very  small  currents  involved.  For  the  garnet  films  on  more 
massive  substrates  there  are  the  additional  questions  of  the 
contribution  of  the  substrate  and  the  interfacial  layer  to 
the  total  conductivity.  When  all  of  this  is  considered  it 


appears  unlikely  that  electrical  measurements  will  find 
much  use  in  the  characterization  of  magnetic  bubble  materials. 

Other  electrical  measurements  which  have  been  carried 
out  in  bubble  materials  include  photoconductivity  in  GdIG 
(Ref. 276)/  the  generation  of  bubbles  by  resistive  heating 
(Ref. 277)  and  the  detection  and  study  of  a  current  controlled 
negative  resistance  in  Si  doped  YIG  (Ref .278  and279 ) .  These 
last  measurements  also  indicate  the  Si  impurity  provides  an 
electrically  active  level  with  an  activation  energy  of 


0. 3eV. 


19 .  SUMMARY 


In  this  report  we  have  presented  the  results  of  a 
study  of  the  physical  properties  of  magnetic  bubble  materials 
Discussions  of  the  known  methods  for  measuring  these 
properties  have  been  given,  and  recommendations,  wherever 
possible,  as  to  which  measurement  techniques  seem  "best" 
at  this  time  have  been  made.  In  order  to  recommend  a 
"best"  technique  among  all  those  which  have  been  cited 
in  the  literature  for  measuring  a  given  material  parameter, 
one  really  needs  to  have  had  experience  with  all  the 
techniques.  Lacking  that,  we  have  considered  such  things 
as  the  simplicity,  accuracy,  and  economy  of  the  various 
methods  along  with  our  experience  with  some  of  the 
characterization  techniques  in  making  recommendations  as  to 
which  methods  are  "best".  Not  all  of  the  physical  properties 
of  bubble  materials  discussed  are  felt  to  be  worthy  of 
routine  characterization.  The  report  has  been  written 
from  the  general  viewpoint  of  someone  just  beginning  to 
set  up  a  bubble  garnet  materials  characterization  facility. 

In  those  laboratories  where  a  strong  capability  already 
exists  in  a  certain  important  measurement  technique,  such 
as  ferrimagnetic  resonance  for  example,  the  "best" 
technique  for  measuring,  say,  anisotropy  field,  at  that 
laboratory,  may  not  be  the  technique  which  we  recommend  here. 
The  results  of  this  study  are  summarized  briefly  in  table  III 
The  reader  is  directed  to  the  appropriate  section  of  this 
report  for  a  more  detailed  discussion. 


TABLE  III 

SUMMARY  OP  RECOMMENDED  CHARACTERIZATION  TECHNIQUES 


Material 

Property 

Frequency 

of 

Measurement* 

"Best" 

Measurement 

Technique  ^ 

Substrate 

Composition 

S 

Electron  Microprobe 

Film 

Composition 

S 

Electron  Microprobe 

Substrate 

Lattice  Constant 

S 

X-Ray  Diffraction  from 

High  Index  Planes 

Film-Substrate 
Lattice  Mismatch 

S 

X-Ray  Diffraction  from 

High  Index  Planes 

Substrate 

Defects 

R 

Etching  and  Microscopic 
Examinations 

Film 

Defects 

R 

Magneto- Optical  Scanning 
in  an  Oscillating 

Bias  Field 

Film 

Thickness 

R 

Interference  of 

Reflected  Light 

Thickness 

Variations 

R 

Photograph  in 

Monochromatic  Light 

Characteristic 

Length 

R 

Strip  Domain  Period 
plus  Thickness 

Magnetization 

R 

Characteristic  Length 
plus  Collapse  Field 

Wall 

Energy 

R 

Derive  from 

Characteristic  Length 
and  4irM 

s 

Wall 

Dynamics 

R 

Strip  Domain  Step- 
Field  Response 

Anisotropy 

S 

AC  Magneto-Optical 

Response  in 

Crossed  Fields 

Coercivity 

i** 

o 

Lock-in  Detection  of 

Partial  Hysteresis 

Loop 

TABLE  III  (Cont.) 


SUMMARY  OF  RECOMMEN DED  CHARACTERIZATION  TECHNIQUES. 


Material 

Property 

Frequency 

of 

Measurement* 

"Best" 

Measurement 

Technique  ^ 

Magnetostriction 

Coefficients 

I 

Strain  Gauge  Technique 
on  Bulk  Samples 

Neel  &  Compensation 
Temperatures 

S 

Magneto-Optical : 

Wall  Susceptibility 
vs .  Temp . 

Temperature 

Variations 

S 

Themoelectric 

Microscope  Stage 

Thermal 

Conductivity 

N 

Use  Literature  Value 
for  YIG 

Index  of 

Refraction 

I 

Brewster  Angle  in 

Monochromatic  Light 

Faraday 

Rotation 

N 

Use  Literature 

Values 

Electrical 

Conductivity 

N 

See  Discussion  in 

Section  18 

*  R:  Routine  Measurement 

S:  Measurement  on  Selected  Samples 

I :  Infrequent  Measurement 

N :  No  Measurement  Should  be  Made 

*t*  • 

See  Discussion  of  Section  19 
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